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955 Abstract—An approximate yet standard procedure, the perturbation method, is used to 
solve analytically the equation for the steafiy state fractionation of multicomponent and complex 
mixtures in non-ideal cascades. It is shown that the formulae are easily adaptable to numerical 
evaluation. The perturbation method is particularly useful when the calculations for a system 
which is close to ideal must be made. The perturbation method may be considered as a non- 
ideality superimposed on the ideal solution. The ideal solution may be thought of as the first 
approximation to the solution while the perturbation is the second approximation. Higher 
order approximations are difficult to carry out. 


Résumé — l’our résoudre analytiquement |'équation relative a la distillation continue de mélanges 
& composants multiples et complexes, les auteurs utilisent la méthode de perturbation, méthode 
classique, bien qu’approchée. 

La méthode de perturbation s'adapte parfaitement aux évaluations numériques : elle est 
particuli¢rement utile quand on étudie un systéme assimilable & un systéme idéal. Elle peut 
étre considérée comme une correction d'état réel sur !hypotheése de la solution idéale. La solution 
idéale correspond & une approximation du premier ordre et la perturbation a une approximation 
du second ordre. Iles difficile d'utiliser des approximations d’ordre plus élevé. 


1. INTRODUCTION be taken into account. If this is done, though 


In this series of articles only ideal cascades with 
constant relative volatilities and constant molal 
overflows were considered, with the exception 
of Part V, where the separation of multicomponent 
and complex mixtures in a packed column was 
investigated. It is found, however, that in many 
instances the mathematical model employed 
so far is too inadequate to represent the actual 
physical system, and that therefore, the variations 
in the relative volatilities of the various com- 
ponents, the unequal molal overflow from plate 
to plate, and the efficiency of the plates must 


the resulting equations become badly non-linear 
and, as consequence, they cannot be solved 
analytically and have been treated in the past 
only by the laborious trial and error plate-by-plate 
procedure, which becomes especially involved 
as the number of the components in the mixture 
to be separated becomes large. 

It is desired here to present a method by which 
the equations of such “ non-ideal ’’ models could 
be solved analytically. The method is only 
approximate and would be expected to be reliable 
only where the deviations from the norm, or 
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ideal, are relatively small; but, in any case, the 
information obtained from these more realistic 
models should prove of value in the difficult 
task of predicting the behaviour of real systems. 

In this article, only variations in the relative 
volatilities and the flow rates with the plate 
number will be considered. However, the method 
can readily be extended to cases where Murphree 
plate efficiencies must be taken into account [1]. 


2. Tue Ipeat CoLtumNn witn VARIABLE 
RELATIVE VoOLATILITIES AND 
UneeualL MoLat OverFLow. THE 
Basic EQUATIONS FOR THE RECTI- 
FICATION OF DiscRETE MIXTURES*® 

The basic equation for the enriching section 
of such a model has already been derived in Part I 


of this series, and is 


R, q(t) 


1 R,, = p,{i) yi) 1+ R, (1) 


Now, strictly speaking, R, and p,{i) are explicit 
functions not of n, but of the temperature and 
the composition of the liquid mixture in each 
stage. It is proposed, however, to circumvent 
this difficulty in the following manner. First 
of all, the model described earlier, that is, in 
which R, R and p,{t) ptt), is completely 
solved, so that the temperature and_ liquid 
composition on each plate can be estimated. 
From this knowledge then, it is possible to 
construct, approximate to be sure but quite 
often reliable, functions R, and p,(i). It will 
be assumed therefore that R, and p,(i) are 
known functions. 

Equation (1) can now be linearized, as was 
done in Part I, by a change in the dependent 
variable. A new function X,(i) is then defined 


so that 


1 alt) X,,(1) 


and therefore 


(8) 


* Notation is the same as in Parts I-IV. 
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or, symbolically, 
= Hy X,(i) (4) 


where H,, is a linear operator. Equation (4) 
can of course be immediately solved by iteration; 
however, since an analytic solution presents 
many advantages, Equation (4) will be solved 
by an approximate yet standard procedure 
known as the perturbation method. This discus- 
sion will be limited only to the first-order 
perturbation, since, from the practical com- 
putational standpoint, the method rapidly loses 
any advantage if the pertubation is carried out 
much further. 

It is assumed then, that X,(i), and H, are 
both functions of a certain parameter z, which, 
in some way, gives a measure of the magnitude 
of the perturbation, and that the problem has 
already been solved for the unperturbed case, 
i.e. when z = 0. X,(i) and H, are next expanded 
in powers of z, so that 


X,(i) = X,(i) +2 +2 Xi) + — — (5) 
and 
H=H+:H+2H + —-— (5a) 


If both the above expressions converge, they can 
be substituted into Equation (4), so that, after 
equating like powers of z, 


k j-k 
X,(i) = 2H X,{i) (5b) 
k=0 


The method of solving Equation (4) then 
consists of solving the unperturbed problem, 
and, by applying (5b) repeatedly with j = 1, 


2,..., of obtaining the functions X,(i), X,(i), 
ete. In this article however, only the first two 
terms in Equation (5) and (5a) will be taken 
into account. 
Let then, with z = 1, 


Pal?) = pli) + 


R, = R + R, (6) 
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It is supposed, moreover, that the perturbation 


terms p,(i), R, and z,(i) are small compared to 
the unperturbed terms. Now since, 


+ pli) j) 
j j 
it follows from Equation (5b), with 7 = 1, that 
= HX ali) + R x a(t) 
+ R, pli) X,(i) + pag (i) 2 27) 
j 


+ pli) X,( 7) 
J 


or 
= A X,(i) + F,(i) + rls, (7) 
where 
X,(i) = X,(i) + X,(¥) (7a) 
and 
fa = xj) Xj) (7b) 
j 


Of course, F,(i) is a known function, because 


X,(i), is the solution to the unperturbed problem 
ond is given by Equation (3-6) of Part I with 


R= _ R and p (i) = pli). 


if rl) is assumed a priori. 

Equation (7) is now a non-homogenious partial 
integrodifference equation which can readily 
be solved. It was in Part I that the eigenfunctions 


Similarly, f, is known 


of the operator H, that is the functions 0,(7), were 
defined by Equations (3-7) and (3-8) again of 


course with = z,(i), pli) = and R = R. 
Therefore, if 


= Cyn) (RA,)"-? v,{i) (8) 


and 


F,{i) = 2 Q(n) v,(t) (8a) 


it follows, from the orthogonality of the eigen- 
functions v,(i) and Equation (7), that 


+1) = C,{n) + + Rf, (9) 
(RR A,)” 
from which it is found that 
n—1 
= + cay (10) 


Finally, since from Equation (3-3) of Part I 


+ p(t) (11) 
R 


X¥,(i) = 


It must be kept in mind of course, that the 


function 2,(¢) and the characteristic numbers 
A, are known, once the unperturbed problem 
has been completely solved. So, the only un- 


known in the solution is the function 2,(#), 
which enters into f,, and which would have to 
be determined from the feed plate matching 
equation. This however can easily be accom- 
plished by a method of successive approximations. 

A very useful asymptotic expression can be 


obtained if 2,(i), for some i, is small. Then, 
compare Equations (8-21)-(3-24), if 
w(i) = 2, hj) (13) 


pili) pj) 

and 

PA) 0) < w (i) 

w(t) 
then 
a, — pli) = 

and 


_ 
pli) 


(% - pli)? 
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Therefore, if z,(i) is sufficiently small 


Xi) = (RA (14) 

Naturally, entirely analogous formulae and 
asymptotic expressions can be obtained for the 
stripping section, but the details will be omitted 
here. 

The advantage then of this analytic solution 
of Equation (4) over the iterative one, lies in the 
fact that the calculations can be performed with 
ease, and the feed plate matching accomplished 
without much difficulty. Also, since the eigen- 
functions v,(i) and the characteristic numbers 
A, are identical with those of the unperturbed 
model, and are therefore independent of the 


function z,(i), the solution is not very sensitive 


to the particular form of z,(i). It must be kept 
in mind however, that the perturbation method 
is only approximate, and, for quantitative 
results, is reliable only when the perturbations 
are small. 


3. THe Extension tro CompLex MIXTURES 


The results of the preceding chapter can now 
be extended to complex mixtures, that is to 
mixtures with an indefinite number of com- 
ponents. Only continuous mixtures will be 
considered here, because the generalization of 
the results to other cases can be carried out 
without much difficulty. 

It was shown in Part I, that, for complex 
mixtures, the basic equations for an ideal column 
contain a distribution function in place of the 
conventional mole fraction; moreover, by 
definition, the distribution function of a contin- 
uous mixture is everywhere continuous, and in 
addition it will be supposed in this chapter 
that said distribution function is everywhere 
positive. 

A question arises however concerning the 
choice of the component designating variable. 
It will be recalled that in the previous articles 
the reciprocal relative volatility p was chosen 
as the component designating symbol. Com- 
plications arise however if, as will be supposed 
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here, the reciprocal relative volatility p varies 
from plate to plate in the cascade. It will be 
assumed again therefore, that the more simplified 
model, that is with constant molal overflow 
and constant relative volatilities, has already 
been solved, so that fairly reliable estimates 
exist as to the composition of the mixture and 
its temperature on each plate. With this informa- 
tion then one could estimate the functional 
dependence of the relative volatilities and the 
flows on the plate number. 

Throughout this chapter then two independent 
variables will be used. The one will be the plate 
designating symbol n and the other will be the 


“component designating symbol” p, where p 
will be the same as the variable p used previously 
in Part I in connection with the “ unperturbed ” 
model. Otherwise the problem to be solved and 
the method of solution parallel those in the 
previous section. 

The basic equation then, for an ideal enriching 
column with variable relative volatilities and 
unequal molal overflow, is 


= _PalP) YAP) + (15) 
pal) ag 
0 


The above can now readily be linearized to 
Xn+i(P) = Ry ilP) 


+ PasslP) | dq (18) 
0 


where the new dependent variable X,(p) is 
defined by 


to P)X P) 
= 17 
Pal P) = ese (17) 
Xn-1 dq 


0 
Again let 
PAP) =P + Pal P) 
R,=R+R, 


= 2p) + 
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and R, = R where it is supposes that the pertur- 


bations rp), R, and 2p) are small in com- 
parison with the unperturbed terms. When 
the above are substituted in Equation (16) 
and only the first order perturbation is taken 
into account, there results that 


X,.:(p) = Rp X,(p) 


ro(q) Xn(q) dg 
+ Rp, 4 X,(P) 


+ R, p X Ap) 


+ Pn4i(P) [ X,,(q) dq 


+ | x n(q) dg 
0 


Finally, Equation 18 can be written as 


=HX(p) + (19) 
where H is the linear operator for the unperturbed 


model — see Equation (4-7) of Part I - and 


j 4 = | 0 X,,(q) dq (20) 


v0 


Again the function F,(p) is known, since it 
contains the solution to the unperturbed pro- 


blem ; and f, is also known if z9(p) is assumed. 

Now to solve Equation (19), use is made of the 
integral transform introduced in Section 4 of 
Part I. It will be recalled, that if f(g) is an 
integrable function of g, then, according to 
Equation (4-18), 


T = p | dg +2 (ay (21) 
0 


where 


A—q 
0 


If then 
T [F,(q)] = 
and 
T (X,(g)] = C,(A) 


Then it is found from Equations (19) and (4-20), 
that 


+ Bf 
(RA)" 


Cy4i(A) = C,(A) 4 


and therefore 


n—l 
R(A)? 


j=1 


C,(A) (23) 


where 


(24) 
So, once C,(A) is calculated, the application of 
the inversion formulae (4°21) and (4-24) yields 


that 


pp COA) (RAM 


A—p 
Cal) 


X,,(p) = 


w( p) 
(p) 
where 


d 


C(A) = - (26) 


[w(A)]? + [A 


Equation (25) then is the approximate solution 
to Equation (16). Since the function w(p) is 
known from the solution of the unperturbed 
problem, the only calculations which have to be 
performed are connected with the numerical 
evaluation of the Cauchy principal value of an 
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improper integral. However, a method for carry- 
ing out such calculations has already been pre- 
sented in Part III. 

It can be said, in conclusion therefore, that 
once the equations for the ideal model, that is 
the ideal cascade with equal molal overflow and 
constant relative volatilities, have been solved, 


then it is possible to extend the results, by the 
method outlined in this article, to any slightly 
non-ideal model. The additional computational 
work needed is shown to be comparatively light. 
The perturbation method, presented here, can, 
moreover, readily be applied to the fractionation 
of any multicomponent and complex mixture. 
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Absorption of carbon dioxide by carbonate solutions in a disc column 


G. H. Roper 
School of Chemical Engineering, N.S.W. University of Technology, N.S.W., Australia* 


Abstract-—Results from the study of the various factors which affect the rate of absorption 
of carbon dioxide by carbonate solutions in a Stephens-Morris disc column are presented. The 
characteristics of the absorption carried out in the disc column were found to be similar to those 
of the absorption in packed towers in respect of the following variables : total concentration of 
the solution, percentage conversion to bicarbonates, rate of liquor flow, and temperature. At 
40%, conversion to bicarbonate, the absorption coefficients for a solution 1-00 total sodium 
normality are given by 
Kg = 182 x (¢ — 23) 


Absorption coefficients for potassium carbonate solutions are approximately 36% higher than 
those for sodium carbonate solutions under the same conditions. 


Résumé—L’auteur présente les résultats d'une étude sur les différents facteurs qui affectent 
la vitesse d’absorption du CO, par des solutions de carbonate dans une colonne a disques 
Stephens-Morris. Par rapport a la concentration totale de la solution, au taux de conversion en 
bicarbonate, a la vitesse d’écoulement de la liqueur et aux températures, les caractéristiques de 
l'absorption dans cette colonne sont semblables a celles de l'absorption dans des tours & garnissage. 

Pour une solution normale en sodium total et avec 40% de conversion en bicarbonate, les 
coefficients d’absorption sont donnés par : 


Kg = 1-82 x 10 (¢ — 23) 


Pour des solutions de carbonate de potassium les coefficients d’absorption sont dans les mémes 
conditions, approximativement 36% plus élevés que pour des solutions de carbonate de sodium. 


INTRODUCTION penetration theory [7] are not in accord with 


Based on the results for the absorption of carbon experimental results obtained on the dise column 
dioxide in water, chlorine in ferric chloride solu- [9], and moreover require a knowledge of the 


tions, and ammonia in water, the Stephens-Morris kinetics of the liquid phase reaction which is 


disc column has been accepted as having charac- generally unavailable and which is difficult to 
teristics similar to those of a commercial packed ™€#sure. If it could be shown that the disc 


tower for physical absorption [12]. Morris and column has characteristics similar to those of a 
Jackson [4] have developed a method for the packed tower for chemical absorption, it would 
laboratory evaluation of liquid film coefficients Probably be possible to develop a laboratory 
for physical absorption in a packed tower based method for the evaluation of liquid film coeffi- 
on data obtained on the same system in a disc cients for chemical absorption based on data 
column. For absorption accompanied by a tained on the same system in a disc column. 


simultaneous chemical reaction in the liquid 

phase, there is at present no method available EXPERIMENTAL APPARATUS 

for the prediction of the liquid film coefficients in AND PROCEDURE 

a full-scale packed tower. The theoretical The StepHens-Morris disc column has been des- 
treatments based on the film theory [5] or the cribed previously [8, 12]. The major dimensions 
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of the absorption unit used in this investigation 
are : 


Number of discs 81 

Dise diameter (mean) 1-475 cm 
Dise thickness (mean) 0-45 cm 
Mean perimeter for liquid flow 0-122 ft. 
Equivalent diameter for gas flow (dry) 0-052 ft. 
Free space (dry) 89% 


Absorption surface (dry) 0-184 sq. ft. 


The discs are threaded edgeways on a 2-mm glass 
road and successive discs are maintained at right 
angles. 

Because of the slow rate of absorption of carbon 
dioxide in the carbonate solutions, the concen- 
trations of the incoming and exit liquor are 
negligibly different for a column of normal 
length. The amount of carbon dioxide absorbed 
cannot be calculated from gas and liquid analyses 
taken at each end of the column, but may be 
calculated from the change with time of the con- 
centration in the graduated 2-litre cylinder used 
as a liquor reservoir. The liquor was pumped 
from the reservoir to a constant head vessel, and 
entered the top of the column after passing through 
a flowmeter (see Fig. 1). Liquid flowed down the 
discs, was withdrawn at the base through an 
adjustable lute vessel, and was returned to the 
reservoir. The liquor in the reservoir was kept 
well mixed by recirculation through the by-pass 
at approximately 330 |b. per hr. The inlet liquor 
temperature was adjusted by means of variable 
electrical heater in the liquor flow line near the 
top of the column. A series of runs consisted of 
taking seven sets of measurements at equal time 
intervals. 

The carbon dioxide was absorbed in sodium 
carbonate solutions (0-31 to 2-16 normal) or 
potassium carbonate solutions (0-86 to 1-41 
normal) at liquor rates from 34 to 415 |b. (hr.) (ft.) 
at temperatures from 76°F to 142°F and at per- 
centages converted to bicarbonate from 2 to 87. 
The amount of carbonate in the liquor was 
determined by titrating a sample from the 
reservoir with standard hydrochloric acid solution 
to phenophthalein, and the total normality by 
titrating to methyl orange. 
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C0; 
BY-PASS 
THERMOMETER 
pump | GRADUATED 2-litre 
CYLINDER 


Fic. 1. Outline of the experimental apparatus. 


The absorption coefficients are expressed in 
terms of the gas phase conditions as given by 


Ko = N,/A- Ap (1) 


where N, is the rate of absorption of carbon 
dioxide, Ib. mol. per hr., A is the surface area of 
the dry discs, sq. ft. and Ap is the driving force, 
atm. The amount of carbon dioxide absorbed in 
one time interval was calculated from the analyses 
of the liquor in the reservoir at the beginning and 
end of the time interval and the total volume of 
liquor in the system, allowance being made for 
the liquor in the lines and other sections of the 
apparatus. The total amount absorbed, M, from 
the beginning of the series is given by the sum 
of the amounts absorbed in each time interval. 
The runs, h, were numbered from — 8 to + 8, the 
central point being O. A fourth-order polynomial 
was fitted to the experimental relation between 
M and h by the method of least squares. The 
constants of the polynomial were evaluated from 
the solutions derived by Avaktan [1]. The slope 
of the curve, dM/dh, at each integral value of 
h, was calculated from the constants by the 
Douglass-Avakian method of numerical differ- 
entiation reported by SHerwoop and Reep [11]. 


ag 
‘ 
» 
4 
VOL 
. 
= 
1955 
3 
> 
= 
td 
q 


Absorption of carbon dioxide by carbonate solutions in a disc column 


Table 1. 
Absorption of carbon dioxide by solutions of sodium carbonate. 


Temp. Ap 
°F 


0-971 

0-971 | 0-327 
0-969 | 0-260 
0-969 | 0-230 
0-964 | 0-202 
0-962 | 0-183 
0-958 | 0-172 


0-970 | 0-146 
0-139 | 0-965 | 0-137 
0-272 | 0-963 | 0-127 
0-397 | 0-959 | 0-118 
0-515 | 0-954 | 0-109 
0-623 | 0-951 | 0-098 
0-714 | 0-948 | 0-089 


0 | 0-956 | 0-190 
0-180 | 0-945 | 0-161 
0-204 | 0-944 | 0-139 
0-472 | 0-940 | 0-122 
0-581 | 0-940 | O-114 
0-740 | 0-937 0-107 
0-752 0-938 0-102 


0-951 0-168 
0-949 | 0-141 
0-045 | 0-125 
0-939 | 0-120 
0-936 | 0-118 
0-934 | 0-116 
0-930 | 0-113 


0-956 | 0-029 
0-029 | 0-954 | 0-027 
0-055 | 0-952 0-026 
0-083 | 0-952 | 0-025 
0-105 | 0-950 | 0-024 
0-129 | 0-949 | 0-023 
0-151 0-948 | 0-022 


0-965 0-234 
0-959 0-202 
0-950 | 0-173 
0-943 0-157 
0-932 0-151 
0-915 0-148 
0-900 0-145 


a 
: 2 i Run Total | | | | | Kg relative 
dM |10° Kg| to 40% 
Series | No. Time concn. | r 
i | dh bicarb. 
No. h min N | | 
| 
1 » 0 1-47 | 341 | 2 |0 3-63 2:79 
—2 45 | «680 | | 1-88 
—-1 | | | 346006 | 81 |: (0-696 | 195 | 1-50 
| 185 | $42 | | 0-962 | 72 | 
1 | 180 | 848 | in | 152 | 41-17 
2 225 | 3460 88 3 | | | 1-38 1-06 
3 270 | | 350 | 4 | 152 | 1-30 1-00 
2 o | 085 | 60-3 79 13 | 109 | 188 
509 | 81 | 20 | 103 | 1-26 
-1 | | 59-6 82 26 0-96 1-17 
he 3: | 0-89 1-09 
Oo | 135 | | 504 | 8 | 8&8 | 
1 | 180 | 59-7 86 | 39 | 0-83 | 1-01 
55 2 | 225 | 60-2 88 44 O75 | 0-91 
3 | 220 60-4 | 89 48 0-68 | 0-83 
3 -3 | | 095 | 181 85 is | 1-44 | 1-57 
-2 | 4 131 | 25 | 1-23 1-34 
1 | 130 % | 32 1-07 | 4 d 
| o | 135 | 31 | 92 | 39 | 0-96 | 
1 | 180 | | 182 | 4 0-68 
coat | 2 | 225 | 133 90 50 | 083 | on : 
ane: 3 | 270 | 132 87 56 | 0-79 | 0-86 
| 
ee 4 —s | o | 188 | 133 78 15 0 128 | 1 “ 
-2 | | 132 so | 19 | | | | 
| 90 131 | 23 | | 096 | 1-409 
| 135 | 87 27 | 0-304 | 0-93 
1 180 130 88 30 | 0-509 | 0-91 
| 2 | 225 | 130 89 84 | (0-625 | 0-90 1-02 
ae 3 | 270 | | 130 9 | #0 | 0740 | 0-88 | 1-00 
| | « 
5 | | 0 | 0-31 | 336 76 | | 0 1 
15 338 | 78 | 45 | | 062 | 
| go | | | | 4 | 0-60 | 
0 45 34-2 79 50 | | 
| 60 844, go | 53 | | oe | 
2 | 75 | 345 | 81 0-52 | 
| 3 90 | 346 | 81 | 61 0-50 | 
| | 76 | 
6 -s | o | 148 | 35 | 73 | | 0 
| | 396 76 | #48 0-202 | | 
9 | | 398 | 78 | 48 | 0-394 1-32 | 
| | 135 3909 | #86 | 55 | 0-556 
1 | 180 | 61 | 0-720 
2 | 225 | 398 82 68 | 0-883 
3 | 20 | 398 | 82 73 «| «108 | 
a | | | | | | 
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Table 1 (cont.) 
Absorption of carbon diowide by solutions of sodium carbonate. 


| Kg relative 
Series No. Time | conen. r Temp. % | 109M) Ap dM | 10° Ke to 40%, 


3 0-180 
; —2 30 | 187 132 | 2 0-162 0-830 0-156 205 | 1-26 
| @ 139 133 35 0312 0821 | 01386 1-79 1-10 
so | 139 134 0 0-455 O811 | 0-121 1-63 1-00 
1 | 120 138 134 43 0-558 0-810 0-109 1-43 0-88 
2 10 | 138 135 47 0-663 0-800 0-102 1-38 | 0-85 
3 | 180 138 135 51 0-762 0-792 | | 14 0-85 


o | 161 304 127 9 0 0-860 0-428 3-61 1-49 
| =§ 45 394 127 18 0405 | 0856 0390 3-30 1-36 
} —2 | 90 304 127 28 0-776 0-852 0-347 | 2-95 1-22 
o | 185 374 128 | 36 | 110 | 0842 | 0808 | 2-65 1-09 
1 | 180 374 129 42 1-39 0-826 | 0269 | 2-36 0-98 4 
2 | 225 374 130 51 1-65 0-811 | 0-236 | 209 0-86 1955 
3 | 270 374 129 53 1-86 0-802 0-206 1-86 0-77 


9 -3 0 143 | 603 14 13 0 0-002 0-207 1-66 1-16 
45 61-4 116 16 0-198 O-804 0-200 1-62 1-13 
—1 | 90 61-4 116 22 0-394 0-892 0-193 1-57 1-10 

Oo | 135 | 61-4 116 |) 29 0-584 0-889 0-187 1-52 1-06 

1 180 61-4 116 34 0-776 0-885 0-181 1-48 1-08 

2 225 61-4 116 38 0-950 0-882 0-176 1-45 101 

3 270 61-4 117 42 1-12 O-874 0-170 1-41 0-99 


0 214 337 115 15 0 0-895 0-276 2°51 117 


—2 40 338 113 24 0-262 0-898 0-257 2-33 1-09 
} —1 | 80 | 338 113 32 0-515 0-896 0-243 2-21 1-03 
i 0 120 339 113 0 0-752 0-889 0-234 2-15 1-00 
3 | 160 338 113 47 0-982 0-886 0-227 2-09 0-98 
} 2 200 337 112 54 1-21 0-876 0-225 2-09 0-98 
; 3 240 | $87 113 62 1-43 0-852 0-220 2-10 0-98 


11 — 3 0 1-87 189 135 20 0 0-834 0-306 2-39 1-03 


45 129 135 25 0-301 0-832 0-301 2:36 101 

me 1 9 | 190 135 81 0-602 0-825 0-298 2°36 101 

; o | 135 190 185 36 0-892 | 0818 O293 | 234 1-00 

1 180 190 135 41 118 | O810 | | 233 1-00 

; | 2 | 225 190 134 46 1-47 0-803 0-285 2-32 1-00 

; | 8 | 270 190 136 51 175 | 0700 | O279 | 2-30 0-99 
3 0824 0479 | B15 

60 400 185 | 30 0444 | 0818 | O420 279 144 

—1 | 120 399 135 38 =| 0833 | 0810 | o373 | 251 | 1-02 

o | 180 | 398 136 45 | 120 | 0790 | 0-343 2-36 0-97 

| 1 | 260 3y9 185 50 154 | 0-787 | 0-827 | 2-26 0-93 

| 2 | 300 390 135 54 | 1-86 | 0-764 | O318 | 229 0-94 

| 3 | 360 | 248 | | 292 0-96 


2 ward. atm dh bicarb 16 
| | | | — 
‘ 
8 | 
10 -8 | 
er 
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Absorption of carbon dioxide by carbonate solutions in a dise column 


Table 1 (cont.) 
Absorption of carbon dioxide by solutions of sodium carbonate. 


| | | | | | 
Run Total | Kg relative 
Series No. Time | conen., T Temp. % 10° M Ap dM |10°Kg| 00% 
No. | h | min. | N | | "F. | bicarb. | atm. dh bicarb. 
| | | | 
3 | -8 0 | 216 482 | 187 | 25 0 | 0-806 | 0-878 | 204 1-51 
2 15 | 48-1 | 136 | 34 | 083 | 0-804 | 9-293 | 1-58 1-17 
a 150 | 48-0 136 42 0-59 | 0-788 | 0-239 | 1-32 0-98 
o | 225 | | 47-9 137 50 | O81 | 0761 | 0-222 1-27 0-94 
1 | 300 | 479 | 137 57 | 106 | 0-732 | 0-214 | 1-27 0-94 
2 375 | 480 | 137 | 65 | 1:26 | 0-682 | 0-206 | 1-81 0-97 


1-37. | 0606 | 0-187 | 1-34 0-99 


0 192 | 372 | 185 | 21 0 | 9-833 | 0-360 | 3-52 1-82 

| 0 135 | 28 | 0336 | 0818 | 0310 | 3-09 1-16 

YL. | =f 80 | 370 135 | 34 | 0602 | 0-810 | 0-278 | 2-75 1-03 
, 0 120 | | 369 185 | 43 | 0-868 | 0-795 | 0-258 | 265 1-00 
155 1 160 | 369 134 | 52 1-14 | O776 | 0-244 | 2-56 0-96 
| 2 |, 200 | | 369 135 60 138 | o743 | 0282 | 255 0-96 

: 240 | | 269 135 | 64 1-49 0-707 | 0-225 | 259 0-97 


15 -3 0 200 =| 39-0 136 16 0 0-820 | 0-360 2:39 1-68 
| | 39-4 136 22 0-348 | 0-815 | 0-290 | 1-93 1-36 

—1 | 120 39-8 135 | 27 0-500 | O817 | O252 | 1-68 1-18 

Oo | 180 41-1 137 31 0-810 0-804 | 0-230 1:55 1-09 

1 | 2 | 40-0 137 35 101 | 0801 | 0-218 1-48 1-04 

300 89-8 187 40 | 121) 0-790 | 0-209 1-44 1-01 


89:8 137 45 1:43 | 0-786 0-202 | 1-40 0-99 


| | | | 
16 | -§ o | 143 | 506 | 1896 | 10 | O 0-821 | 0-283 | 2-50 1-48 
_2 45 | sos | 196 | 18 0-266 | 0-820 | 0255 | 225 1-33 
eS 90 59:1 136 | 26 0-510 | 0-816 | 0-230 2-04 1-21 
0 135 | 58-8 136 33 0-729 | 0-810 | 0-209 | 1:87 
1 180 585 | 136 39 0-914 | 0-805 | 0-189 | 1-70 1-01 
58-2 136 45 1:10 0-794 | 0-175 | 1-60 0-95 
127 | | 0-168 155 | 0-92 


| 
17 —3 0 118 | 67-9 i386 | 5 | 0 0-825 | 0-328 | 288 | 202 
fi @ I | 63-8 | 135 | 18 | 0-80 0-822 | 0-261 | 2-30 1-62 
| -1 90 | | 697 | 185 | 28 | 052 | 0-820 | 0208 | 184 | 1-80 
0 135 | we | 185 | 87 | O72 | 0812 0-172 | 153 | 1-08 
1 180 | | 706 | 185 43 | 085 | 0807 | O146 | 131 | 0-92 
| 70-6 | 135 50 1-01 0-795 | 0-126 | 1-15 | 0-81 
| | 55 0-784 | 0-108 1-00 | 0-70 
| 
| | o | 1-82 | 415 | 134 1 | 0 | 0534 | 461 | 168 
-3 45 | 413 136 | 20 | 0-487 | 0-624 | 0-429 | 377 | 1-87 
| —1 9 | | 411 | 137 30 «| 0-882 | 0814 | 0360 | 3-21 1-17 
135 | 408 138 37 1-21 0-804 | | 2-82 1-03 
180 408 | 135 43 1-52 0-804 | 0-296 | 2-67 0-97 
| 135 50 1-82 0-795 | 0-284 | 2-59 0-94 
| | 0-267 


. 
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Absorption of carbon dioxide by solutions of potassium carbonate. 


| T 
Run | | Total | | Kg relative 
Series No. Time | r | concn. % Temp. | 10° M Ap dM | 0° Kg to 40%, 
No. h min | N bicarb. "F atm. | dh bicarb. 
| 

1 | 397 | 1-41 | 88 | | 0-962 | O389 | 2-93 1-29 
—2 45 307 | 18 83 0-371 | 0962 | 0362 | 273 1-20 
-1 | 395 | 30 | 81 | o719 | 0-962 | | 2-41 1-06 
0 135 | 893 40 82 | 104 | 0959 | 0304 | 2-20 1-01 
1 180 392 | 51 82 0-955 | 0-282 2-14 0-04 
2 225 | 302 | 62 | 82 | 1461 0-950 | 0275 | 208 0-91 
3 0-037 | 0-269 0-01 
0 | | | 0-954 | 0-260 1-95 
—2 45 | 546 | 12 87 | 0232 | 0-951 | | 157 | 1:56 
—1 9 | 528 | 20 87 0-418 | 0-951 | 0-174 1-33 | 1-32 
0 135 50-9 | 27 89 | 0-580 | 0-948 | | 120 | 119 
1 180 | 506 89 | O731 | 0-948 | O142 | 1-08 1-07 
2 225 | 303 42 88 | 0-870 | 0-948 | 0182 | 101 1-00 
8 100 | 0950 | 0-125 0-95 
3 -—8 0 | 48:8 | 0-99 4 136 || 0 0-819 | 0255 | 3-89 1-54 
-2 30 48-9 | 14 | 186 | 0244 | 0-819 | 0-232 | 3-08 1-40 
-1 60 49-0 | 23 | 136 | 0-464 | 0-817 | 0209 | 278 1-26 
0 90 49-1 | 32 | «(136 | 0-684 | 0815 | 0186 | 2-47 1-12 
1 120 40-4 39 136 | 0-870 | 0813 | O162 | 2-17 0-99 
2 150 49-7 46 142 | 101 «| O-779 2-02 0-92 


-3 30 | 401 | 

0 | 401 53 
1 | 120 | 401 68 
| 10 | 401 79 | 
8 | 180 | 401 87 


136 


136 | oO 0-817 | 0410 | 5-46 1-40 
136 | 0-871 | 0815 | 08380 | 4-40 | 1-18 
132 | 0-661 | 0-829 | 0-282 370 | 0-95 
185 | 0940 | O811 | 0256 | 343 | 0-88 
131 1-22 0-811 | 0-234 | B14 0-81 
133. | (1-45 | 0-222 | 8-11 | 0-80 
133 | 1-62 0-722 | 0-159 | 239 | 0-61 


1-28 


0-805 


| 

0-145 
0-128 

| 


The rate of absorption is given by dividing the 
slope by the length of the time interval in hours. 
These lengthy calculations were undertaken to 
eliminate any personal error in measuring the 
slope form a graph of M against h. The partial 
pressure of carbon dioxide over the carbonate- 
bicarbonate solution was determined from 
equations given by Perry [6, p. 679]. The partial 
pressure of water vapour of the solution was 
taken to be the same as that of a sodium carbonate 
solution of the same total normality, and was 
determined from data of Perry [6, p. 172]. The 
partial pressure of carbon dioxide in the gas 


The experimental results for the absorption of 


RESULTS AND Discussion 


phase was calculated by subtracting from the 
total pressure in the column the partial pressure 
of water vapour. The density of the solution was 
calculated from data of Perry [6, p. 181]. 
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carbon dioxide in sodium carbonate solutions are 
given in Table 1, and for absorption in potassium 
carbonate solutions in Table 2. The absorption 
coefficients have been correlated in terms of a 
standard solution 1.00 total normality 40% con- 
verted to bicarbonate. 


Table 2. 
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Absorption of carbon dioxide by carbonate solutions in a disc column 


EFFECT OF 
CONVERSION TO BICARBONATE 


For each series, the readings were taken at almost 
constant values of total normality, liquor rate, 
and temperature. The effect of the conversion to 
bicarbonate was obtained by plotting the experi- 
mental coefficients as a function of the percentage 
bicarbonate in the solution. The coefficients at 


40%, bicarbonate was read from the graph, and 
the relative coefficients were obtained by dividing 
the experimental coefficients by the figure. The 
liquor rate, temperature, and driving force, Ap, 
were plotted on the same graph, and values at 


— COMSTOCK AND DODGE 


RELATIVE COEFFICIENT (Fp) 


70 
% BICARBONATE 


Fic. 2. The effect of the absorption coefficient on the 
conversion to sodium bicarbonate (one third of the points 
are plotted). 


40% bicarbonate noted. The values of the co- 
efficients, liquor rates, and temperatures read 
from these graphs are given in Table 3 and are 


Table 8. Absorption of carbon dioxide in carbonate solutions. 


used in the correlations showing the effect of 
liquor rate, temperature and total normality. The 
coefficients relative to those at 40% bicarbonate 
are shown in Fig. 2 for sodium carbonate solutions, 
and in Fig. 3 for potassium carbonate solutions. 


| | 


COMSTOCK AND DOOGE 


2 


RELATIVE COEFFICIENT (Fg) 


| Lt 
10 
% POTASSIUM BICARBONATE 


Fic. 3. The effect on the absorption coefficient of the con- 
version to potassium bicarbonate. 


Comstock and Dopce [2] have measured the 
rates of absorption of carbon dioxide into 
potassium carbonate solutions over a wide range 
of conversions to bicarbonate. They used a 8in. 
diameter tower packed with glass rings 10 mm 
o.d. x 8mm id. x 10mm length to a height 
of 10-67 ft. The relative coefficient corresponding 
to the arithmetic mean of the per entage bicarb- 
onate in the liquor entering and leaving the tower 
has been calculated from their data and is shown 


Data at 40% conversion to bicarbonate. 


Series Temp. 


10° Kg | °F. 


| 


Data for sodium carbonate solutions 
1-30 
0-82 
0-92 
0-88 
0-65 
1-65 
1-63 
2-42 
1-43 


60-0 


33-6 
395 
139 
374 

61-4 


own 


117 


(b) Data for potassium carbonate solutions 
1 2-28 82 393 
2 101 88 5 


om 
| 
ae © | Total | Series | | Temp. | Total 
No. r | comen.N | No. | | °F. | concn. N 
(a) | | 
1-47 10 214 | 118 839 1-26 
0-85 11 2338 | (185 190 1-87 
, eas | 0-95 12 2-45 185 | 394 1-99 
ai 1-38 13 1-35 136 48-0 2-16 
aoe 0-31 14 2-66 135 369 1-92 
Ss | 1-43 15 1-42 137 39-8 2-00 
oa 1-19 16 1-69 136 58-4 1-43 
ae 1-61 17 1-42 135 70-6 1-18 
pets | 1-49 18 2-75 136 408 1-82 
1-41 8 2-20 136 49-4 0-99 
Pee 0-86 4 3-90 135 401 0-95 
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by the lines on Fig. 2 and Fig. 3. These lines 
represent satisfactorily the data from the disc 
column. 


Errect or Liquor Rate 


The correlation between the absorption co- 
efficients and liquor rate, temperature, and total 
normality was made by the method of successive 
approximations. The final correlation is 


Kg = 1-82 x 10° F,- F,(t — 28) r* (2) 


The effect of liquor rate is shown in Fig. 4. Those 
series taken at or near 85°F and 135°F were 


mak 


| 


| 


| 


Fic. 4. Correlation of the absorption coefficients with 
liquor rate for sodium carbonate solutions corrected to 
1-00 total normality and 40% conversion to bicarbonate. 


corrected for temperature to either 85°F or 135°F 
and for concentration to 1-00 total normality. 
The equations for the lines on Fig. 4 are 


Kg = 0-118 x 10° (3) 
at 85°F, and 
K, = 0-204 x 10° (4) 


at 135°F. A different slope at the two tempera- 
tures would fit the data better, but the difference 
between the slopes is small considering the scatter 
of the data. 

Furnas and BEe.uincer [3] have studied the 
rates of absorption of carbon dioxide in sodium 
carbonate solutions, using a 12in. diameter tower 
packed with gin. and lin. Raschig rings and 
lin. Berl saddles. For physical absorption, the 
liquid film coefficients are proportional to the 
0-54 power of the liquor rate for jin. rings, to 
the 0-76 power for 1 in. rings, to the 0-65 power 


for 1 in. saddles [10], and to the 0-7 power for 
dises [8, 12]. It is possible that the liquid film 
coefficients in a packed tower can be related to 
those in a dise column by 


L™ (5) 


where n is the slope of the coefficient-liquor rate 
line plotted on log-log paper for the absorption 
when carried out in the dise column, and z is the 
ratio of the slope for physical absorption in a 


6 


ee 

| 

Yein RASCHIG RINGS | | 
DATA OF FURNAS AND BELLINGER 
| 
| | 


l 


Ib/(hrisq ft) 


Fic. 5. Data of Furnas and Bevuincer for j in. Raschig 


rings compared to the predicted slope. 


packed tower to the slope for physical absorption 
in a disc column. On this basis, 2 is 0-54/0-7, or 
0-77, for Zin. rings, 1-09 for 1 in. rings, and 0-98 
for lin. saddles. The corresponding slope for 
the absorption of carbon dioxide in sodium 
carbonate solutions are 0-47 x 0-77, or 0-36, for 
Z# in. rings, 0-51 for 1 in. rings, and 0-44 for 1 in. 


4 


A 


RASCHIG RINGS 
DATA OF FURNAS AND BELLINGER 


ib/(hedisq ft) 


Fic. 6. Data of Furnas and Bevurncer for 1 in. Raschig 
rings compared to the predicted slope. 
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Absorption of carbon dioxide by carbonate solutions in a disc column 


The data of Furnas and BELLINGER 
has been 85°F, 
1-00 total normality, and 40% conversion to 
The data for 
g in. rings are shown in Fig. 5, for 1 in. rings in 


saddles. 
recalculated and corrected to 
bicarbonate using equation (2). 


Fig. 6, and for 1 in. saddles in Fig. 7. The data 
in Fig. 5 scatter badly and no conclusion can be 
drawn, but the data of the other two figures 
show that the predicted slope is in good agreement 
with the experimentally determined slopes. 

4 


DATA OF FURNAS AND BELLINGER 


| | | 


1000 10000 


Ib/(hrX\sq ft) 
Data of Furnas and Bevuincer for 1 in. Berl 
saddles compared to the predicted slope. 


Fic. 7. 


The data of Comsrock and Dopce [2] using 
glass rings cannot be correlated in this manner 
due to the unknown effect on the liquid film 
coefficient for physical absorption of the small 
diameter of their column, and of the large height 


Ox 


= 


ft) 


Fic. 8. Data of Comstock and DopcGe showing the effect 
of liquor rate at various temperatures on the absorption 
coefficient using potassium carbonate solutions. 


of packing without redistribution of the liquor. 
Their data have been recalculated to 1-00 total 
normality and 40% conversion to potassium 
bicarbonate at temperatures of 59°F, 104°F and 
133°F, and is plotted in Fig. 8. Part only of their 
28 runs at temperatures near 59°F is shown, but 
all the runs near 104°F and 133°F are shown. 
Within the accuracy of the data and the same 
slope is applicable at each of the temperatures. 
Errect oF TEMPERATURE 

The experimental coefficients were corrected for 
liquor rate, total normal normality and con- 
The coefficient at 85°F 
was taken as unity, and the coefficients relative 
to 85°F calculated. The data are shown in 
Fig. 10. The points at 85°F and 135°F are taken 
from the lines of Fig. 5, and therefore represent 
a large number of individual coefficients. The 
points from the data of Comstock and Dopce 
[2] were taken from the lines of Fig. 9, and 
represent a large number of individual coefficients. 
The data of Fig. 9 can be correlated by one line, 
showing that the effect of temperature on the 


version to bicarbonate. 


absorption in a dise column is in close agreement 
with the effect of temperature in a packed tower. 
22 


x COMSTOCK AND DODGE 
@ THIS INVESTIGATION 


“a 


RELATIVE COEFFICIENT ——e 


100 120 140 
TEMPERATURE 


Fic. 9. Effect of temperature on the absorption coefficient. 


Errect or Tora. CONCENTRATION 


The effect of total concentration is small, almost 
within experimental error. Comstock and 
Dopce [2] have found that an increase in total 
concentration decreases the coefficient slightly. 
The data of this investigation were divided into 
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sections of low, medium and high concentrations. 
The coefficients corrected for temperature, liquor 
rate, and conversion to bicarbonate, were aver- 
aged in each section. These data are shown in 
Fig. 10. The line representing the results of 
Comstock and Dopce is shown on the same 
figure. Within experimental error, the effect of 
total concentration is the same for both columns. 


COMSTOCK AND DODGE 

~ 

= THIS INVESTIGATION 

~ 
= 

l i 

1-2 6 2-0 24 


TOTAL NORMALITY ——> 
Fic. 10. Effect of total normality on the absorption 
coefficients. 


COMPARISON OF SOLUTIONS OF 
SoDIUM AND PoTassiumM CARBONATES 


In addition to their extensive data for potassium 
carbonate solutions, Comstock and Dopce give 
data obtained using sodium carbonate solutions. 
Compared under similar conditions of tempera- 
ture, total concentration and conversion to 
bicarbonate, their coefficients for potassium 
carbonate solutions are approximately 40% 
higher than their coefficients for sodium carbonate 
solutions. The coefficients for the dise column 
given in Table 3 show that those coefficients for 
potassium carbonate solutions are 36% higher 
than the coefficients for sodium carbonate 
solutions. 


G. H. Roper 
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CONCLUSION 


The dise column has been compared to packed 
towers for the absorption of carbon dioxide by 
carbonate solutions. The disc column shows 
characteristics similar to those of packed towers 
for the absorption in respect to the effect on the 
absorption coefficient of the liquor rate, the total 
normality of the solution, the percentage con- 
version to bicarbonate and temperature. 
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NOTATION 


a = surface area of dry packing, sq. ft. per cu. ft. 
A = surface area of dry discs, sq. ft. 

Fp = factor to allow for the conversion to bicarbonate (i.e. 
the coefficient relative to the coefficient at 40% 
bicarbonate) 

Fc = factor to allow for the effect of total normality (i.e. 
the coefficient relative to the coefficient at 1-00 total 
normality) 

h = the number of a run in a series, taken at equal time 
intervals 
ky, = the liquid film coefficient of absorption, ft. per hr. 

Kg = the overall coefficient of absorption based on gas 
film conditions, Ib. mol./(hr.) (sq. ft.) (atm.) 

Ky = the overall coefficient of absorption based on gas 
film conditions, Ib. /(hr.) (sq. ft.) (atm.) 

L = rate of liquor flow in a packed tower, Ib./(hr.) (sq. 
ft.) 

M = the amount of carbon dioxide absorbed, Ib. mol. 

N = total normality 

N4 = the rate of absorption of carbon dioxide, Ib. mol. 
per hour 

Ap = the driving force forabsorption, atm. 

t = temperature, °F. 
I’ = rate of flow of the liquor, Ib./(hr.) (ft.) 


{1} Avaxtan: Thesis in mathematics, M.I.T., 1933, reported by Saerwoop and Reep, p. 266. 
[2] Comstock, C. W. and Dopnce, B. F.; Industr. Engng. Chem. 1953 29 520. 

{3} Furnas, C. C. and Betiincer, F.; Trans. Amer. Inst. Chem. Engrs. 1938 34 251. 

[4] Morris, G. A. and Jackson, J.; Absorption towers, London, Butterworth, 1951. 

(5) Peaceman, D. W.; Thesis in chem. eng., M.I.T., 195i, 206-31. 

(6) Perry, J. H.; Chemical Engineers’ Handbook, 3rd edn., McGraw-Hill, New York, 1950. 
(7] Perry, R. and Picrorp, R. L.; Industr. Engng. Chem. 1953 45 1247. 


(8) Roper, G. H.; Chem. Eng. Sci. 1953 2 18. 
{9} Roper, G. H.; Chem. Eng. Sci. 1953 2 247. 


{10} Smerwoop, T. K. and Hotioway, F. A. L.; Trans. Amer. Inst. Chem. Engrs. 1940 36 39. 
{11]) Suerwoop, T. K. and Reep, C. E.; Applied Mathematics in Chemical Engineering, 1st edn., McGraw-Hill, 


New York, 1989 287. 
[12] 


Sreruens, E. J. and Morais, G. A.; Chem. Eng. Progress 1951 47 232. 


264 


q 4 
4 
1 
| a 
¥ 
VOL. 
4 
1955 
| 
qj 
: 


Chemical Engineering Science, 1954, Vol. 3, pp. 265 to 268. Pergamon Press Ltd. 


The effect of a surface-active agent on mass transfer in a stirred 
liquid-liquid extractor* 


A. and S. G. Tersesen 


The Technical University of Norway, Chemical Engineering Laboratory, Trondheim, Norway 


Abstract—The effect has been observed of small additions of sodium oleyl-p-anisidinesulphonate 
on the rate of extraction of iodine from aqueous solutions with carbon tetrachloride in an 
agitated extractor. The results have been analysed by the method of Hrxon and Smrru [2]. 
The time required for extracting 90% of the extractable iodine from the aqueous phase increased 
with the concentration of surface-active agent to a maximum about 70% above the original 
value at 3 x 10-5 g/100 ml aqueous solution, and fell on further additions far below the original 
value. The Hixon lines were curved for the larger additions of surface active agent, indicating 
a gradual increase in interfacial area. 


Résumé—Les auteurs ont étudié l’effet produit par l’'addition de petites quantités de sulfonate 
de sodium d’oleyl-p-anisidine, sur la vitesse d’extraction de l'iode & partir de solutions aqueuses 
au moyen de C Cly, dans un extracteur agité. Les résultats ont été analysés par la méthode de 
Hrxon et Smiru. Le temps nécessaire 4 l'extraction de 90% de liode que peut fournir la phase 
aqueuse, croit avec la concentration de l'agent tensio-actif jusqu’é un maximum de 70% environ 
de la valeur primitive de 3 x 10° g pour 100 moles de solution aqueuse. Les courbes de Hrxon 
s‘infléchissent pour des additions plus importantes d’agent tensio-actif, ce qui indique une 


augmentation graduelle de la surface a l'interface. 


INTRODUCTION 
In recent years several investigators have studied 
the effect of surface-active agents on the rate of 
Both 
spray and packed columns have been used, as 
The most important 


mass transfer in liquid-liquid extraction. 


well as agitated systems. 
papers on the subject have been reviewed recently 
by Garner and Hace [3] who have themselves 
found that the rate of extraction of diethylamine 
from toluene by drops of water was reduced to 
45% of the original value by the addition of 
0.015% Teepol. They were able to show that this 
retardation appeared also with droplets which did 
internal circulation. Although the 
not 


not exhibit 


literature is unanimous, a_ considerable 


amount of evidence has been accumulated to 
show that surface-active agents, adsorbed at an 
interface between two liquids, can effect a sub- 
stantial reduction in the rate of mass transfer. 
There is, however, still insufficient data available 
to enable firmly based conclusions to be drawn 


with regard to the nature of this effect. 


° A summary of this paper has appeared elsewhere [1]. 


The present work was carried out as a pre- 
liminary to further investigations. A stirred 
tank extractor was used, and the results were 
analysed by the method of Hrxon and Smiru 
[2]. The system studied by these authors, viz. 
carbon tetrachloride — water — iodine, was chosen 
also for the present work, and sodium oleyl—p- 
anisidinesulphonate was selected as surface-active 
agent. 


EXPERIMENTAL 


The apparatus was in principle similar to that 
used by Hrxon and Samira [2], and is illustrated 
in Fig. 1. 

A cylindrical glass vessel of about 4 litre total 
capacity was equipped with a simple paddle 
stirrer driven at 140 + 1 r.p.m., and arranged 
with its lower edge at the interface as shown. The 
vessel was placed in a water bath which was 
thermostatically controlled at 22-0°C. In all 
experiments the volume of the aqueous phase was 
10 times that of the organic phase. 


+ Present address : Notodden Salpeterfabrikker, Notodden, Norway. 
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Fic. 1. Extraction apparatus. 


Iodine of standard analytical quality was 
dissolved in distilled water to give a concentration 
of about 230 mg/l. 3 litres of this solution was 
added to the apparatus, which was covered with 
paper to reduce loss of iodine by evaporation. 
The stirrer was started, and a sample taken. 
800 ml of refractionated, commercial carbon 
tetrachloride was added carefully through a 
funnel which reached nearly to the bottom of 
the vessel. The extraction proces; was followed 
by photometric determinations of iodine in 
samples withdrawn from the vessel. During the 
runs only the aqueous phase was sampled, but 
at the end the organic phase was also analysed 
as a cheque on the material balance. The volume 
of the samples withdrawn was compensated for 
as described by Hrxon and Smiru. 

With no addition of surface-active agent the 
action of the agitator caused wave-formation at 
the interface, but very few droplets of carbon 
tetrachloride were detached from the main 
bulk. With increasing additions, more free 
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droplets were observed, and with still larger 
additions a froth appeared at the interface. In 
no case was more than a small fraction of the 
organic phase dispersed. 

The surface-active agent was obtained as 
Lissapol LS, a commercial product of the J.C.J. 
Ltd., Dyestuffs Division. This was dissolved in 
warm n-butanol, filtered, crystallized several 
times, and finally washed with ether and dried. 
For the present experiments a stock solution was 
made containing per 100 ml : 

0.1023 g sodium 
0.1656 g sodium sulphate 
0.0225 g sodium chloride 


Stock solutions older than 14 days were discarded. 
In the present paper all concentrations of surface- 
active agent are expressed as g of the active 
component per 100 ml of the aqueous phase. 


RESULTS 


On the basis of the two-film theory of extraction, 
Hixon and Smiru [2] have deduced that a plot 


We 
of log (1 — W = against time should give a 


straight line, provided the interfacial area and 

the partition coefficient remain constant through- 

out arun. The value of Pearce and EveRsoLe 

[4] for the partition coefficient at 25°C, 

= = 99.7 78 iodine per mi carbon tetrachloride 
mg iodine per ml water 

has been used in the present work. This gives 

= 0-893 x W°,p, which enables the 

Hixon curves to be plotted. The results of some 

typical runs are shown in Fig. 2. 

It can be seen that the lines are straight except 
for the larger additions of surface-active agent. 
The deviation from the straight line indicates an 
increase of the interfacial area with time, and 
makes it impossible to characterize a run by a 
single constant as Hrxon and Smiru have done. 
Instead, the time ty, required for extracting 90% 
of the extractable iodine (Wg = 0-9 x Weg) 
was noted as well as 2 tg 4, where is the 
time required to give Wy, = 0-684 x W,,,. For 
a straight line these figures should be the same, 
and the difference can be taken to give a measure 
of the departure from linearity. 


Ca 
Wee 
: 
= 
| 4 
1955 
dss 
4 
i 
; 
7 
7 
i? ig 
7 


Effect of a surface-active agent on mass transfer in a stirred liquid-liquid extractor 


\ = | CONCENTRATION OF SODIUM OLEYL-p-ANISIDINE - 
-SULP T 
© NO ADDITION ~SULPHONATE 


2% 10% 
272 16% 
° 


O08 


40 60 80 
EXTRACTION TIME MIN ——= 


Fic. 2. Typical runs with and without surface active 
agent, plotted according to Hrxon and Sirs [2]. 


EXTRACTION TIME, MIN 


oO a 5 6 7 8 9 
CONCENTRATION OF SODIUM OLEYL- p- ANISIDINESULPHONATE, ¢/ ml AQUEOUS PH 


Fic. 8. Time for extracting 90% of extractable iodine, 
ty and 2 X tgg4 plotted against concentration of surface 
active agent. 
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The results have been summarized in Fig. 3, 
where fy, and 2 X tg, are plotted against the 
concentration of surface-active agent. 

With increasing concentration, the extraction 
times increase to a value about 70% above the 
original at 8 x 10° g/100 ml, and subsequently 
fall well below the original. These results are 
analogous to those obtained by Cuu, TayYLor 
and Levi [5] from measurements of the rate of 
extraction of benzoic acid from benzene into 
water in a packed column. At small concen- 
trations of weak surface-active agents there is a 
retardation of the mass transfer due to the 
molecules adsorbed at the interface, but at 
higher concentrations the tendency of the agent 
to increase the interfacial area predominates. At 
about 3-7 x 10% g/100 ml a sudden decrease in 
the extraction time resulted from a smal! addition 
of surface, active agent. This critical concentration 
was found to vary greatly with the speed of the 
agitator. 

These results differ from those already reported 
in the literature by the extremely small con- 
centrations of surface-active agent found to be 
effective in reducing the rate of mass transfer, 
1 part per 10’ parts of water being sufficient to 
The effective con- 


give a noticeable effect. 
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centration is thus about 500 times smaller than 
that reported for Teepol by Garner and Hate 
[3]. It can be seen from Fig. 3 that t,, and 
2 X tgg.4 are equal for concentrations of surface 
active agent less than about 2 x 10° g/100 ml, 
which means that the theory of Hrxon and Surru 
holds in this range. This result indicates that the 
process responsible for the retardation of mass 
transfer is established more rapidly than the 
extension of the interfacial area caused by the 
larger concentrations of surface-active agent. 


Acknowledgement—The authors wish to thank 
the Royal Norwegian Research Council for a grant 
which enabled this work to be carried out. 


NOTATION 


m = Partition coefficient : 
mg iodine per ml of solvent phase 
mg iodine per ml of aqueous phase 


Wy, = Weight of solute in extract phase (organic) g 


We eq = Weight of solute in extract phase (organic) after 
equilibrium is established g 
Wp° = Weight of solute in raffinate phase (aqueous) at 
the start of the experiment K 
fgg = Time required for extracting 90 %, of the extract- 
able iodine (Wp = 0-9 WE eq) min 

teg.4 = Time required to give Wp = 0-684 Wp, min 


[2] Hrxon, A. W. and Smirn, M. L; Ind. Eng. Chem. 1949 41 973. 
[3] Garner, F. H. and Hauer, A. R.; Chem. Eng. Sci. 1953 2 157. 


[4] Pearce, J. N. and Eversore, W. G.; J. 
Cuv, J. C., Tayvor, C. C. 


Phys. Chem. 1924 28 245. 
and Levy, D. J.; Ind. Eng. Chem. 1950 42 1157. 
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A simple method of measuring gaseous diffusion coefficients 
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(Received 28 April 1955) 


Abstract—In this paper a new and simple method of measuring gaseous diffusion coefficients 
is described. The method was developed so as to enable coefficients to be readily obtained with 
sufficient accuracy for design purposes. The method consists essentially in measuring the amount 
of equi-molar counter-diffusion occurring in a known time along a tube of known bore and length 
connecting two reservoirs initially containing gas mixtures of different concentrations. 

Values of diffusion coefficients obtained by this method are compared with previous measure- 
ments for two of the gases, and satisfactory agreement is found. Coefficients are given for the 
diffusion of CO,, NHg, SO, Cl, and Brg in air at 20°C and 760 mm Hg pressure. 


Résumé—L auteur décrit une méthode simple et nouvelle pour mesurer rapidement et avec 
une précision suffisante, les coefficients de diffusion gazeuse. La méthode comporte essentielle- 
ment la mesure de la diffusion 4 contre-courant d'un nombre égal de molécules, dans un tube 
(section et longueur connues) nécessitant deux réservoirs contenant au départ, des mélanges 
gazeux de concentrations différentes. L’auteur trouve ses résultats en accord avec ceux des 
auteurs précédents. I] donne les coefficients pour la diffusion de CO,, NHg, SOg, Cl, et Brg 
dans l’air & 20° et sous 760 mm. 


When using this apparatus for measuring the 
diffusion coefficient of, for instance, a certain gas 
in air, one of the flasks is initially filled with a 
mixture of air and the gas, the stopcock in the 
diffusion tube being closed as also are the connec- 
tions to the absorber. The other flask is filled with 
air to approximately the same pressure. If the 
gas is lighter than air, the mixture is introduced 
into the upper flask, if heavier the lower flask is 
filled with the gas mixture. The density difference 
thus opposes convective transfer between the 
flasks when the stopcock in the diffusion tube 
is opened. 


APPARATUS AND PROCEDURE 


A knowledge of the gaseous diffusion coefficient 
for a specified gas mixture is frequently necessary 
for the design of absorption and stripping towers. 
Reference to published tables of physical pro- 
perties too often shows the desired diffusion 
coefficient to be missing. The method of measure- 
ment to be described was developed so as to 
enable gaseous diffusion coefficients to be readily 
obtained with sufficient accuracy for design 
purposes. 

The apparatus used, sketched in Fig. 1, is of 
very simple construction, consisting of two 
flasks of about 2 litres capacity connected by a 
glass tube of about 0-7 cm bore and 24cm long 
having a stopcock of the same bore in the centre. rapid 
The flasks are set up with connecting tube constant temperature throughout the mass. 
vertical in a large (150 litres) unheated water In order to eliminate any bulk flow through the 
bath. Other connections to the flasks, shown diffusion tube which would occur when the 
in the sketch, permit filling and purging either stopcock is opened if the gas pressures in the 
flasks and also connection of the flasks to an _ two flasks were unequal, both flasks are connected 
absorber vessel. The latter consists of a horizontal to the absorber for a period prior to opening the 
tube partially filled with an absorbent liquor or stopcock in the diffusion tube. The absorber 
granulated absorbent depending on the gaseous acts as a trap for the gas contained in the air/gas 
diffusion coefficient being measured. mixture and ensures that only air can be 


After filling, the whole apparatus shown in the 
sketch is left in the water bath for several hours, 
stirring ensuring the attainment of a 
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transferred from one flask to the other during 
the final pressure equalising operation. Both flasks 
are then isolated from the absorber and the 
stopcock in the diffusion tube is opened. 


Equimolar counter diffusion of the gas and air 
along the diffusion tube is then allowed to proceed 
for from one to five days, after which the stopcock 
is closed and the contents of each flask are 
analysed. 


Sketch of diffusion apparatus. 


Fig. 1. 


INTERPRETATION OF RESULTS 


The process taking place while the stopcock is 
open is essentially that of equi-molar counter 
diffusion of the two gas species along a tube of 
known length and cross-section resulting from the 
concentration difference between the two ends 
of the tube. 


Thus if 1 = Effective length of diffusion tube 
== true length + end corrections. 


a = Cross-sectional area of diffusion 
tube. 


Q = Quantity of either species trans- 
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the volume of pure gas at the 
total pressure and temperature in 
the apparatus. 

D = Diffusion coefficient under prevail- 
ing conditions of temperature 
(T° K) and pressure (p) in the 
apparatus. 

Initial concentration difference of 
one species between the two flasks, 
concentration being expressed on 
a volumetric ratio basis. 


AC, = Concentration difference at time 
(t) of this species between the two 
flasks. 


Then ignoring the transfer associated with 


setting up the initial linear concentration 
gradient along the tube, at any instant :— 

dQ Da 

1 


Now if v, and v, are the volumetric capacities 
of the two flasks then :— 


(2 
dQ 


Substituting for dt from (1) in (2) gives :— 


(AC,) = Da (* + ~) dt 


Integrating from AC,= AC, to AC, = AC, 
and t = 0 tot =¢ then :— 
l 1 AC 
aii. 1). AG, 


AC, 
_ (xc) 
~ at (AC, — AC,) 


then :— D 


The length (l) is the effective length of the 
diffusion tube and comprises the true length 


ferred in time (t), expressed as plus two end corrections which allow for the 
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resistance to transfer associated with spreading 
of the diffusing species from the ends of the 
tube throughout the volume of the flasks. Each 
end resistance is virtually equivalent to the 
resistance to transfer from one side of a circular 
flat disc into a surrounding infinite medium. 
This case has been solved mathematically [1], 
the resistance being equivalent to a uniform 
layer of fluid on the disc having a thickness of 
md/8 where d is the disc diameter. Hence :— 


1 = True length of diffusion tube + md/4. 


For the apparatus used, the end correction 
amounted to about 2% of the effective length 
of the diffusion tube. 

In deriving the above expression for calculating 
the diffusion coefficient, the state 
transfer necessary to set up the initial linear 
concentration gradient which follows the discon- 
tinuous gradient existing before opening the 
stopcock, has been ignored. The unsteady 
state transfer associated with the changing 
concentration profile in the tube during the run 
has similarly been ignored. The volume of gas 
contained in the diffusion tube is relatively 
quite small compared with the total transfer 
during a run, and with runs of over one day’s 
duration it can be shown that the maximum 
error resulting from neglecting the unsteady 
state effects is less than 4% in the value of the 
diffusion coefficient. 


unsteady 


A simple method of measuring gaseous diffusion coefficients 


Table I. Measurements of Diffusion Coefficient of Dilute CO,-Air. 


Some CoEFFICIENTS MEASURED 


In order to effect a comparison between the 
values of diffusion coefficients obtained using 
the simple apparatus and those measured 
previously by other investigators using more 
conventional methods, a series of measurements 
was made of the diffusion coefficient of dilute 
CO, in air at about normal temperature and 
pressure. The results of six independent measure- 
ments with the simple apparatus are given in 
Table I, the values of the diffusion coefficient 
being reduced to 20°C and 760 mm Hg pressure 
by assuming that the coefficient is inversely 
proportional to the pressure and proportional 
to the 7/4 power of the absolute temperature. 

The principal dimensions of the diffusion 
tube and flasks employed for the above measure- 
ments are : 


True length of diffusion tube 23-4 em 
Mean cross section of diffusion tube 0-410 cm? 
Volumetric capacity of upper flask 2300 cm* 


Volumetric capacity of lower flask 2278 cm*® 


The above six runs give a mean value of 
0-165 cm*/sec. with a mean deviation of 0-008 
em*/sec. Recent work by Boyp et al, [2]. using 
a conventional double cell apparatus, but employ- 
ing an interferometer for following the unsteady 
state diffusion, gave a value of 0-160 cm?/sec. 
for the diffusion coefficient of dilute CO* in N, 
at 20°C and 760 mm Hg. Previously BoarpMaAN 


t hours 89-66 63-75 


| 
Pp mm Hg 755 760 | 


64-67 66-0 | 
755 765 | 


55 : 
665 | 111-5 
745 | 765 
| 203 203 | 293 291 | 292 291 
we 
| 0-1132 0-1108 | 0-1121 0-1136 | 0-1074 0-1035 
| 
AC, — | 00788 0-0632 | 0-06-46 0-0638 0-0589 0-0884 
| | 
Q em® | 39-4 | 54-5 54-4 56-9 55-5 74:5 
ek D em®/sec. | 0-170 0-164 | 0-164 0-159 0-170 0-167 
en 
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and WiLp (1987) using a double cell apparatus 
obtained a value of 0-163 cm*/sec. for the same 
system. 

A series of twelve measurements using runs 
of from 16 to 70 hours duration have also been 
made to determine the diffusion coefficient of 
dilute NH, in air. These gave a mean value for 
the coefficient at 20°C and 760mm Hg of 
0-227 cm,/sec. with a mean deviation of 
0-008 cm,/sec. This may be compared with 
0-225 cm,/sec. obtained by Wintercetst [4). 

Other coefficients measured with the apparatus 
are :— 


Dilute SO,- Air D,,° = 0-122 cm* /sec.: mean 
deviation 0-005 cm? sec. 

Dilute Br,-Air = 0-091 cm*/sec.; mean 
deviation 0-004 cm? ‘sec. 

Dilute Cl,-Air = 0-124 c¢m*/sec.; mean 
deviation 0-002 cm? /sec. 


REFERENCES 
{1] McApams, Heat Transfer McGraw-Hill, 2nd Edition, Table II. 
[2] Boyp, Srey, et al. ; J. Chem. Phys. 1951 19 548. 
[3] BoarpMan and Witp ; Proc. Roy. Soc. 1937 A162 511. 
[4] Wryrercetst ; Ann. Phys. 1930 4 3 323. 
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No previous measurements of these three coeffi- 
cients have been traced. 

It appears from the above figures that this 
type of apparatus can be used for obtaining 
reproducible values for gaseous diffusion coeffi- 
cients. Values are in satisfactory agreement 
with those measured by other means where 
this has been done. 


NOMENCLATURE 


a = Cross-sectional area of diffusion tube. em? 
D = Gaseous diffusion coefficient cm? /sec. 
i Effective length of diffusion tube em 
p = Total gas pressure in apparatus mm Hg. 
Q@ = Volume (at pressure p and temperature 7) 
of either species transferred during run em® 
{ = Duration of diffusion run sec, 
T = Gas temperature in apparatus °K 


AC, = Initial concentration difference of one 
species between the two flasks 

AC, = Concentration difference at end of run of 
one species between the two flasks 
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Extraction of acetic acid from water 


4. Ternary-vapour-liquid equilibrium data 


F. H. Garner, S. R. M. and C. J. Pearce 


Devartment of Chemical Engineering, The University, Birmingham, 15 


Abstract—This paper presents ternary vapour-liquid equilibrium data for the system ethyl 


acetate-acetic acid-water. 


The data correlates satisfactorily, and linear interpolation can be 


used to predict the ternary vapour liquid equilibria from the binary data when an appropriate 


correction factor is used. 


Résumé—Les auteurs presentent des résultats relatifs a l'équilibre liquide-vapeur du systéme 


ternaire acétate d’éthyle-acide acétique-eau. 


La corrélation des données est satisfaisante et il est possible, en tenant compte d'un facteur 
de correction convenable, de prévoir par une interpolation linéaire, l'équilibre ternaire a partir 


des données binaires. 


INTRODUCTION 

Tue system ethyl acetate-acetic acid-water is of 
great commercial interest but there has been 
little information published on ternary vapour 
liquid equilibria ; accurate data on this subject 
are required both in the analysis of the operating 
conditions of azeotropic columns and for the 
design of plants of higher efficiencies. 


EXPERIMENTAL 


The reagents used in the determinations were of 
the same purity as those described in previous 
papers, Part 1, 2 and 3. 

When both the vapour condensate and the 
liquid phase in equilibrium were single phase, the 
equilibrium still of ELi1s [24] was used. When the 
liquid and condensate phases were partially 
immiscible, the still described by SELVARATNAM 
[34] was used. 

A binary mixture of acetic acid and water was 
made up and was diluted with ethyl acetate after 
each determination. 

Experimental results were thus determined 
along a line of constant molal ratio of acetic acid 
to water. After allowing three hours for equili- 
brium to be established, the liquid and vapour 
condensate were sampled into ice-cooled bottles. 
When a single phase was present, samples were 
then analysed immediately for water and acetic 


acid content by titration with Fischer reagent 
and barium hydroxide respectively. When 
ethyl acetate was present in small quantities the 
ethyl acetate was estimated by saponification 
with potassium hydroxide solution. 

If either the vapour or liquid condensate was 
two phase, a known weight of alcohol was added 
to produce miscibility. This method was prefer- 
able to separating the two liquid phases and 
analysing them separately, as the samples were 
small and often difficult to separate. 


RESULTS 


The vapour-liquid equilibrium data are listed in 
Table 14 and are represented in Figs. 14 to 19. The 
activity coefficients were calculated according to 
the equation : 


(1) 
wh -re y, = activity coefficient of component 1 
nm = total pressure 
mole fraction of component 1 in the 
vapour 
- mole fraction of component 1 in the 
liquid 
= vapour pressure of component 1 at 
the temperature of the determination. 
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The above equation applies strictly only to 
solutions in equilibrium with an ideal vapour 
mixture. The correction factors are uncertain 
for the acetic acid vapours due to the association 
of the acetic acid vapour, and are smal! for ethy! 
acetate and water, so no correction factors have 
been used. 

The experimental results in Figs. 14, 15 and 16 
are plotted by the method suggested by Hanns 
and Norman [86]. For a mixture of components 
A, Band C the procedure used is as follows. The 


ratio —— 


AaB in the vapour is plotted against the 


corresponding ratio 7uB™ the liquid, using 


the molal concentration of the third component 


100 B 
C, as a parameter. Similar ratios of Bue and 


100C . 
quem the liquid are plotted against ratios 
100 B 100 C 
B+ 
meters A and B respectively. Figs. 17, 18 and 19 
represent the log activity coefficients of 4, B and 
C plotted against the mole percentage of A, B and 
100B 1004 1004 
B+CA+CA+B 
respectively as parameters. Lines of constant 
parameters were then interpolated between the 
experimental points as suggested by CoLBURN 
and SCHOENBORN [5]. 


in the vapour using the para- 


used 


C in the liquid with 


CORRELATION OF RESULTS 


Methods for prediction of ternary data from 
binary data and for correlation of scanty ternary 
data, have been suggested by numerous authors 
[18], [88], [89], [40], and [41]. These methods of 
correlation and prediction may be classified into 
four main groups as follows : 

(1) Linear interpolation. 

(2) Volume fraction interpolation. 

(3) Calculation of ternary activity coefficients 
by thermodynamically consistent two, three and 
four suffix equations, involving a number of con- 
stants calculated from binary data and, in some 
cases, ternary constants calculated from experi- 
mentally determined ternary data. 
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(4) The use of equations expressing ratios of 
the log activity coefficients (e.g., log y,/ye, 
log y2/7g, etc.) in terms of the mole percentages 
of each component in the liquid, and constants 
obtained from equations expressing the same 
ratios from the binary systems. This involves 
the use of ternary constants in certain cases. 

In applying these methods to the irregular 
system ethyl acetate-acetic acid-water, we are 
concerned with the prediction of vapour liquid 
equilibria in both the miscible and immiscible 
ternary regions. 

In the miscible ternary region the following 
difficulties arise. 

(a) The correlation of the two binary systems in 
which acetic acid is one of the components 
becomes complicated by the association and 
inter-association effects which are possible 
with this component. In a ternary system it 
might be expected that there would be an 
even more irregular behaviour. As no really 
successful method has been found for the 
correlation of acetic acid-water the ternary 
results have been calculated using a molecular 
weight of 60 for acetic acid. 

(b) Fig. 9, Part 3 Chem. Eng. Sci. 1954 3 48, 
shows that the binary system acetic acid- 
water cannot be correlated by the use of 
van Laar or Margules equations. 

(c) Fig. 11, Part 8, Chem. Eng. Sci. 1954 3 48, 
shows that in the immiscible region of the 
system ethyl acetate-water, the van Laar 
equation does not hold. The equation is of 


K 
the form log y = —. In a ternary mixture 
z 


the mutual solubility of ethyl acetate and 
water is increased by the presence of acetic 
acid. The problem therefore arises as to 
whether or not in this increased miscible 
region the van Laar equation should be used 
for predicting ternary data. 


From a consideration of the above points it is 
unlikely that equations of the type previously 
mentioned under groups (8) and (4) will be 
accurate, since if the ternary equations are to be 
applicable, the binary systems must obey binary 
equations of a similar form. It would thus appear 
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that these methods may not be applicable to 
systems exhibiting partial miscibility or inter- 
association. 

On the other hand, linear interpolation methods 
have the great advantage that it is not necessary 
for the binary data to be correlated by equations 
of similar form. Such methods are basically 
limited to certain symmetrical systems but it will 
be shown that by the introduction of a correction 
factor the method is applicable to irregular 
solutions. 

From Figs. 4, 5 and 6 and from the data of 
Pratr [42] and work in our laboratories [43] it 
can be seen that linear interpolation does give a 
reasonable approximation to the experimental 
results. The accuracy of linear interpolation to 
irregular systems can be improved by the intro- 
duction of an appropriate correction factor. 

The proposed equation is of the form 


og = z, V1) (logis ¥1) (4) 


where log y, = log activity coefficient of com- 
ponent 1 in the ternary system at concentration 2, 
in the liquid. 

= log activity coefficient of component 
1 in the binary system 12 at some concentration 
of z, of component 1. 

logis ¥; = log activity coefficient of component 
1 in the binary system 18 at some concentration 
2, of component 1. 

Under certain conditions the final term of the 
above equation may become very small in which 
case linear interpolations between the binary 
systems will give a good approximation to actual 
behaviour. 

Using a reduced form of the two suffix Margules 
equation and rearranging the equation, the 
following expression was obtained : 


2s 
+ 22, (Ag, 
+ [Ag — Ay, — Agg — 22, 

(Ag, — + 225 

(Agg — Ags) — C(1—22,)) (5) 


logy, = [ + 22, (Ag, — Aj,)} (1—2,)? 


Ajs)| (1+2,)* 


+ 


It will be noted that in this equation the terms 
[Ais + 2a, (Ag, — (1 —a,)? and [Ais + 
2x, (Ay, — Aj 3)](1 — 2,)* are the values of log y, 
in the binary systems at a concentration z, of 
component 1. These two terms in the above 
equation are therefore equivalent to the value 
given by linear interpolation between two binaries 
which obey the three suffix Margules binary 
equation. A similar analogy can be made for the 
2 and 4 suffix Margules equations. 

The assumption is now made that the final 
term is approximately of the same value even if 
the binary systems do not obey the Margules 
equation. The following equation is therefore 
proposed for irregular solutions such as ethyl 
acetate-acetic acid-water. 


log y; [logis “i Je, 
8 


7 [logys Yi Je, 


— Ajyg— Agg—22, (Ag, — 
+227,( Ass Ags) — Cc —2z,)] (6) 


+ 


The constant C has to be obtained from ternary 
data, and in the system ethy] acetate-acetic acid- 
water, C = 0, has been used in equation (6). 


COMPARISON OF METHODS 


In Table 15 the values of activity coefficients 
predicted from binary data, using three different 
methods of prediction are compared with experi- 
mental ternary activity coefficient data (column 2) 

Column 8 gives the results for linear interpolation 


using the activity coefficient plot, log y = - for 


ethyl acetate-water in Fig. 11. If the van Laar 
correlation is used for ethyl acetate-water the 
discrepancy between the predicted and experi- 
mental results is greater. 

Column 4 gives the predicted activity coefficient 
for equation (6) using the binary activity coefficient 


plot log y = : in Fig. 11 and the activity coeffi- 


cient plot for acetic acid-water in Fig. 9. 
Column 5 gives the predicted activity coefficients 


based on the three suffix Margules equation. 
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Table 15. 


Water-acetic acid-ethyl acetate. 


Predicted from 
3-suffia Margules 
0 


Predicted from 
eq. (6) C 0 


Predicted by 
linear interpolation 


Experimental 


results 


Liquid 
compositions eqns. C = 


YA YE 


2-02 
1145 
1-64 
1-65 
4-65 
0-955 


10 


0-78 
0-782 
0-56 
0-94 
0-86 
0-922 


0-71 


= 


a 


A 


o obtain 


It will be seen that the values given by equation ternary liquid mixture is immiscible. new 


(6) are much nearer the experimental values than type of still is being designed in order t 


The results in further experimental data in this region. 


either of the other two methods. 


column 4 calculated from equation (6) are shown It is suggested however that ternary vapour 


to be within 5-0%, of the experimental values. compositions in the immiscible ternary region can 


It is interesting to note that in the increased — be predicted from equation (6) and the use of the 
miscible region (see 4 above under Correlation of - :' , 
equation log y for the immiscible binary 


Results) the closest agreement between predicted 


system. The vapour composition can be pre- 


and experimental results is obtained by equation 


dicted either from liquid compositions based on 


(6) and the log y plot for ethyl acetate-water. 


This 


an overall composition in the immiscible region 


unexpected conclusion is being further or from ternary compositions on the solubility 


investigated. 


PREDICTION OF IMMISCIBLE VAPOUR 
LiqeuIp EQUILIBRIA 


Only a limited amount of experimental work has 


curve. 
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Abstract—The properties of an ideal solution are carefully reviewed, and their application 
to actual solutions of gases and liquids is considered. A number of general equations relating 
the properties of solutions of gases and liquids, both ideal and non-ideal, are developed. 


Résumé—L auteur fait avec soin la revue des fonctions caractéristiques d'une solution idéale 
et de leur application aux solutions réelles de gaz et de liquides. 


Il expose des équations générales reliant les propriétés de gaz et de liquides idéaux ou non 


idéaux. 


INTRODUCTION 


THE concept of an ideal solution originated many 
years ago, and one might expect that the implica- 
tion of the term would by now be quite clear to 
all who have studied thermodynamics. Unfor- 
tunately, this does not seem to be the case, at 
least among engineers. Thus, a careful re-examina- 
tion of the ideal solution and its properties would 
seem to be justified. At the same time it is 
appropriate to consider the deviations of actual 
solutions from ideality. Engineers make use of 
such quantities as the volume change of mixing, 
the heat of mixing, the entropy change of mixing, 
etc. These quantities are related, though not in 
a simple way, and in theory at least they may 
all be calculated from experimental measurements 
of the volume change of mixing. The required 
relationships are derived below. This treatment 
will be limited to solutions of liquids and gases in 
which no dissociation or chemical reaction occurs. 

Three conclusions are usually drawn from the 
statement that a solution is ideal : 


(a) Vi=V, V;=V;, ete, and 


therefore AV = 0, as may be seen by com- 


paring the equations : 
V (1) 
V + A4V (2) 


For gases this is equivalent to AMaGat’s law. 


H, =H, 4,=H, H,;=H,, ete., and 
therefore AH = 0, as may be seen by com- 


paring the equations : 


i’ 


H = 2,1. +%H, +... 
H = 2H, + %H, + 


+ 2,H,, (3) 
+2,H, + AH (4) 


This is equivalent to the statement that the 
integral heat of mixing is zero. 


(c) The Lewis and Randall rule is valid, ice. : 

I, = (5) 
Eqs. (1, 3) and (5) apply either to gaseous or 
liquid solutions, where all components are stable 
when pure in the same physical state as the 
solution at the temperature and pressure of the 
Eqs. (2) and (4) are obviously com- 
pletely general since AV and AH are defined so 
as to make them correct. 

It is not my purpose to dispute these equations, 
for they certainly are correct for an ideal solution 
as originally defined. The original definition was 
clearly stated by Lewis and RanpDa tu [2], who 
described an idea] solution as one in which the 
fugacity of each component is proportional to its 
mole fraction at all concentrations. This imme- 
diately leads to eq. (5), the Lewis and Randall 
rule. In addition, they stated that this pro- 
portionality, i.e., eq. (5), must hold at all 
temperatures and pressures. It is this part of 


solution. 
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the definition which is sometimes overlooked. 
Indeed, it is not to be found in any of the standard 
textbooks on chemical engineering thermo- 
dynamics. The complete definition is usually 
repeated in textbooks on thermodynamics for 
chemists, but the necessity of the latter require- 
ment is seldom made clear. More important, 
some explanation is required as to what is meant 
by all temperatures and pressures, particularly 
with regard to liquid solutions. 

An incomplete understanding of the ideal 
solution has led some engineers, at least, to 
erroneous concepts. Notable among these is the 
notion that a solution is necessarily ideal, or 
nearly so, if at a given temperature and pressure 
its volume change on mixing is zero or small at 
all concentrations. That this is not so will be 
shown presently. 


2. DERIVATION OF 
GENERAL EQUATIONS 


Let us now proceed to the development of general 
equations relating the quantities being discussed, 
starting with the relation between fugacity and 
volume. These derivations are for a closed system 
of constant mass and constant over-all or total 
composition. The system is not limited to a 
single phase, but may consist of both liquid and 
vapour. As a matter of convenience we will base 
all equations on one mole of material, so that the 
total amount of any component i will always be 
given by its mole fraction in the entire system, z,. 
By the first law of thermodynamics for a closed 

system, 
dE = dQ —dW (6) 


For any reversible process in this closed system 
in which the only force is that of fluid pressure, 


dQ = TdS (7) 
and 

dW = pdV (8) 
Therefore, 

dE = TdS — pdV (9) 


This equation relates properties only, and is 
therefore divorced from any consideration of the 
process causing the change even though a 
reversible process was assumed in arriving at it. 


By definition, 
and 


dF = dE + pdV + Vdp — TdS — SAT 


Substitution for dE by eq. (9), gives, 


dF = Vdp — SdT (11) 
At constant temperature, 
dF = Vdp (12) 


This familiar derivation has been repeated here 
to emphasize the fact that no limitations need be 
imposed concerning the number of components 
or phases. Thus F, E, V, and S refer to the total 
system, in this case one mole, regardless of its 
complexity. A consequence of eq. (11) is that 
p and T are always independent variables for F 
for any system whatever. 

Eq. (12) may be written: (dF /dp)y = V. 
Differentiating with respect to N, the number 
of moles of component é : 


dp) 
dN; 
Where the symbol z-- indicates constant total 

composition. This may also be written : 
oF 
dp 


By the definition of a partial molal quantity this 
becomes : 


(oF, 


or, at constant temperature and composition, 
dF, = V, dp (13) 


F,, the partial molal free energy of component i 
in solution, is also known as the chemical potential 
and is often given the symbol ,,. 


Again no assumption has been necessary as to 
the number or nature of the phases. Integration 
of this equation from p*, the zero-pressure state, 
to p at constant temperature and composition, 
gives : 


| (14) 
pe 
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A further equality may be added to this equation 
through the definition of fugacity for a com- 
ponent in solution. 


Thus : 
F, =| Pap RT -\n(J,/z,p*) (15) 


In a similar fashion eq. (12) may be integrated 
for pure component i at the same temperature 
and pressure, and combined with the defining 
equations for fugacity of a pure material, to give : 


p 

F,—F; =| = RT - \n(f;,/p*) (16) 
p* 

Subtraction of eq. (16) from eq. (15) gives : 


P 

F,— F,—(Fj — Fj) -| (V,— Vi) dp 
Since [1], F; — Fj; = RT -In(z,), the above 


equation becomes : 


— F, — RT -\n(z,) = — V,)dp 
p® 


= RT -\n(f;/z,f,) (17) 


It must be understood here that in the above 
equations corresponding values of F, and F,, 
f, and f,, and V;, and V; are always taken at the 
same temperature and pressure in their actual 
states at these conditions even though the pure 
component may be in a different physical state 
than the mixture. 

Eqs. (15, 16) and (17) are usually represented as 
applying only to single-phase systems ; however, 
such a limitation is not at all necessary, and 
indeed is incorrect. These equations are derived 
from the first law of thermodynamics together 
with the definitions of free energy and fugacity. 
The only assumptions required are that the 
system be closed and that the only force be that 
of fluid pressure. As seen from the derivations, 
no assumptions are necessary concerning the 
number or physical state of the phases either for 
a solution or for the pure components. Thus, 


integration may be carried through the two-phase 
region, and indeed must be, if these equations are 
applied to the liquid phase. 

In connection with eq. (17) it should be pointed 
out that the group, f;/z,f, is commonly called 
an activity coefficient and that its deviation from 
unity may usually be regarded as an indication 
of the deviation of the solution from ideality. 
Such an interpretation is correct only when f; 
represents the fugacity of component i in a 
single-phase solution and f, represents the fugacity 
of pure 7 in the same physical state as the solution 
at the temperature and pressure of the solution. 
For the latter, fictitious states are sometimes 
employed. In the integration of eq. (17), only 
actual states are used, and as a result the term, 
f:/%f is properly regarded as an activity 
coefficient only for the single-phase regions at 
conditions such that the pure components can 
exist in stable form in the same state as the 
solution at the given temperature and pressure. 

In eqs. (15, 16) and (17) the integrations are 
carried out at constant temperature, and constant 
over-all or total composition. The symbol, z,, is 
used here to represent these over-all or total 
compositions so as to distinguish them in the 
two-phase region from the liquid and vapour 
compositions which are usually represented by 
x, and y,; respectively. Of course, for a single 
phase z, becomes identical with z; or y,, but in 
the two-phase region they are clearly different. 
It must be understood that while total com- 
position is held constant for the integrations of 
eqs. (15) and (17), the phase compositions must 
vary in the two-phase region in accordance with 
the requirements of the phase rule since the 
temperature is also constant. In this connection 
a few general comments about the phase rule 
will not be amiss. 

The only variables considered in this regard 
are temperature, pressure, and the individual 
phase concentrations. Total volume, total mass, 
total mass of any constituent, etc. are not phase- 
rule variables. When we speak of fixing the state 
of a system, we do not imply anything with res- 
pect to its size. Thus the phase rule may be used 
to predict what changes can be made in an open 
system without changing its state. But we are 
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here everywhere considering a closed system in 
which changes occur under restraints quite differ- 
ent from those in an open system. If we fix the 
total amounts of each constituent (as has been 
done by fixing the total composition and restricting 
the system to one mole of material), the only 
remaining variables are temperature and pressure, 
and fixing them determines all properties of the 
closed system provided that temperature and 
pressure can both vary, and vary independently. 
Thus we may write for such a system, G = ¢(T, p) 
where G represents any property such as enthaply, 
entropy, or free energy. 
not be varied or 


If temperature and 


can not vary 


pressure can 
independently, then this relation fails. 
an exception to this, for if G is the free energy 
(F), then this relation is always valid because 
F = 4(T, p) under all circumstances as shown 
by eq. (11). 


In the above equations, V, and F, are partial 


There is 


molal quantities, defined in the usual manner. 
These quantities are commonly thought of in 
connection with single phases, but they certainly 
have unique values for two-phase (or even 
multi-phase) systems at equilibrium, and the 
integration of eq. (17) through or into the two- 
phase region assumes the use of these quantities 
for the total system in the two-phase region. This 
point does not seem to have been discussed before, 
and therefore deserves elaboration. We have long 
been conditioned to think of fugacity and the 
partial molal quantities as applied only to single 
phases, but this is an unnecessary restriction. 
Consider the physical path analogous to the path 
of integration. We start with a vapour mixture 
of definite composition and definite amount at a 
pressure approaching zero, and compress at 
constant temperature and constant composition 
and of material. For 
pressure the system will have a definite set of 
properties even after the dew point has been 
reached. Indeed, for any pressure between the 
dew point and the bubble point, there will be a 
definite division of the material between the two 
phases, and each phase will have a definite com- 
position : 2,, 2, etc. in the liquid and y,, y, ete. 
in the vapour. Thus for each pressure all properties 
of the two-phase mixture are fixed, and volume, 


constant amount each 


enthalpy, entropy, free energy, and the partial 
molal properties all have definite values. In the 
two-phase region the values of f; and of F, or 
#,; must be the same for the liquid phase as for 
the vapour phase. Nothing about eq. (17) is 
contrary to this requirement. The equations 
developed in this paper do not in any way conflict 
with equations in common use, but rather 
supplement them. 

Let us turn our attention now to the relation 
between the volume change which occurs on 
mixing and the integral heat of mixing. It has 
sometimes been assumed that these quantities 
are directly related, and that if one is zero or 
small, the other must also be zero or small. While 
this is by no means true, the two quantities are 
related, though not in any simple fashion. The 
equation relating them may be developed most 
directly as follows: (The derivation is for a 
binary mixture, but it may obviously be general- 
ized for any number of components.) From 
eq. (4): 

AH 


=H 


-2,H, —%H, (18) 


For a mixture [2], 
H = H* + | (V—T(dV/dT] dp 


p® 


(19) 


Again the integration is carried out at constant 
temperature and total composition, and again the 
equation applies not only to the vapour, but also 
throughout the two-phase and liquid regions. The 
only requirement for the application of this 
equation to a closed system is that H be a 
function of both temperature and pressure. This 
can be true only when temperature and pressure 
are independently variable. In the case of a single 
phase, this is always true. For a two-phase 
system (the maximum we are considering here), 
temperature and pressure are not independently 
variable for a single component or for an azeotrope 
in multi-component systems. But for mixtures 
other than azeotropes, this equation is valid under 
all conditions in single- and two-phase systems. 
For azeotropes a latent-heat term would have to 
be added if the integration were carried through 
the two-phase region. 
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While eq. (19) may be used in general to 
calculate enthalpies of both liquid and vapour 
mixtures, it can not be used as it stands to 
calculate the enthalpies of pure components in 
the liquid region. For a pure component, con- 
densation at constant temperature also occurs at 
constant pressure, and the pressure and tem- 
perature are not therefore independent variables 
in the two-phase region. For this case a latent- 
heat term must be added to eq. (19) to take 
into account the enthalpy change caused by 
condensation at temperature and 
pressure. Thus we may write for the pure com- 
ponents, a and b: 


constant 


= + [iv 
ge 


p 


H, = 4 | [V, dp —L, (20b) 


-T(IV,/dT)| dp —L, (20a) 


L, and L, are the mola! latent heats of pure a and 
Eqs. 
(20a) and (20b) are written specifically for the 


pure 6 at the temperature of the mixture. 


liquid phase, but they are valid for the vapour 
phase if the latent-heat terms are omitted. 

Substitution of eqs. (19), (20a), and (20b) in 
eq. (18) gives : 


AH = H* — +2,L, —2%H, 


+ 2,7 + 2,T (=)| dp 

V —2z,V, —2V, = AV, and 

—(V/IT) + 2,(3V,/IT) +%(dV,/2T) = 
-(JAV/DT). 


But by eq. (2), 


Therefore, 


+ — — + + 
But by eq. (4), 


However, AH*, the integral heat of mixing at p*, 
approaches zero for a value of p* approaching 
zero. 

Therefore, 


p 


AH - [AV — T(JAV/dT)] dp +2,L, + z,L, 


Generalizing, 


T (AV dp + {zL,) (21) 


As it stands, eq. (21) gives the relationship 
between the integral heat of mixing and the 
volume change of mixing for liquid mixtures. If 
the latent-heat terms are omitted, it applies to 
gas mixtures. Clearly, the heat of mixing depends 
not only on the volume change of mixing at the 
temperature and pressure of mixing but on the 
volume changes at all lower pressures as well. The 
integral in eq. (21) is evaluated at constant 
temperature and total composition. It should be 
emphasized that for a liquid mixture this integra- 
tion must be carried through the two-phase region 
and that AV always represents the difference 
between the mixture volume and the volumes of 
the pure components in their actual states at the 
temperature and pressure of the mixture even 
though these states may be different. Of course, 
AV will ordinarily be finite for actual gases at a 
pressure approaching zero even though the ideal- 
gas law is valid at this pressure. 

We can now complete the picture by deriving 
several more relationships between the A’s of 
mixing for solutions of both gases and liquids. By 
eq. (17) (constant temperature and composition): 


F, — F, — RT - in(z,) (V, — V,) dp 


But 


Therefi re 


H,—TS,;—H,+TS,—RT -\n(z) = V;) dp 


p® 
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Multiplying through by z; and summing over all 


components : 


RT &(z,\nz,) = | [2 (2,V,) — (z,V,)) dp 
By the definitions of AH, AS, and AV, this 
reduces to: 
AH — TAS — RT 2 = [ (avydp 
But at constant temperature, AH — TAS = AF: 


Therefore 


AF = + RT &(z, In z,) (22) 
pe 


By the definition of the excess free energy of 
mixing, i.e., the free energy change over and 
above that for mixing ideal gases, 


AF,, = AF — RT 2 (z, Inz,) 


Therefore 
AF., = [avrap 23) 


Since AS = (AH — AF)/T, then by eqs. (21) 
and (22) : 


v 
AS = — (QAV/dT) dp — R 
” + 2(z,L,)/T (24) 
By the definition of the excess entropy of mixing, 


i.e., the entropy change over and above that for 
mixing ideal gases, 


AS,, = AS + R (z, Inz,) 


Therefore 


= @avprydp + 2(2,L,)/T (25) 
pe 


Eqs. (24) and (25) as written apply to liquid 
solutions. For a gaseous mixture the latent-heat 
term must be omitted. The integrals in these 
equations are, of course, evaluated at constant 
temperature and total composition. 
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8. APPLICATION OF THE GENERAL 
E@avuATIONS TO IDEAL SOLUTIONS 


The equations developed above are general 
relations applying alike to ideal and non-ideal 
solutions. For the case of the ideal solution, 
however, they may be greatly simplified. 

Let us apply eq. (17) to the case of an ideal 
solution as defined by Lewis and Ranpa.u. For 
such a solution, f; = 2,f,; therefore In ( f,) 
= 0, and eq. (17) reduces to : 


F.—F RT In (z,) (26) 


This equation is a general relation for a given 
solution only if true over a range of temperatures 
and pressures. This accounts for the requirement 
that f, = z,f, for an ideal solution at all tem- 
peratures and pressures, for we are about to 
differentiate eq. (26) first with respect to pressure 
and then with respect to temperature, and such 
differentiations would not be admissible were this 
not a general equation. For an ideal solution 
as originally defined, it is completely general, as 
are the equations which result when it is differ- 
entiated. The requirement that the Lewis and 
Randall rule holds at all temperatures and 
pressures obviously restricts consideration to a 
hypothetical solution in which no phase changes 
occur. 

Differentiation of eq. (26) with respect to 
pressure gives : 


(OF ,/dp) — (F,/dIp) = 0 


Therefore {1}: V, = V, and AV =0 for an 
ideal solution. 
Dividing eq. (26) by T and differentiating with 


respect to temperature gives : 


(F/T) d(F,/T) 
oT oT 


= 0 


Therefore —H,/T* +H,/T* =0. Thus 
H, = H, and AH = 0 for an ideal solution. The 
same conclusion is reached by substituting 
AV = 0 into eq. (21) and noting that this must 
hold at all temperatures and pressures. Since 
this can be true only if no phase change is 
involved, the latent-heat term must be omitted. 


In the same way, eqs. (23) and (25) may be used 
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to show that AF,, and AS,, are zero for ideal 
solutions. 

It should be clear from the above discusson 
that an ideal solution may be defined equally 
well by stating that it has any one of the three 
attributes listed at the beginning of this paper, 
provided that this attribute is stated to hold at 
all pressures, temperatures, and concentrations. 
This is all very satisfying from the mathematical 
point of view, but the engineer is interested in 
using these relations, as approximations at least, 
for real solutions, and no real solution can be 
ideal at all temperatures and pressures. The 
practical problem then is one of determing when 
or over what range of conditions these equations, 
true for ideal solutions, may be used as good 
approximations for real solutions. 


OF GENERAL EQUATIONS 
SOLUTIONS 


4. APPLICATION 
To Non-IpDEAL 


Consider now a few general cases of non-ideal 
behaviour. They will be represented by diagrams 
of (V, — V,) vs. pressure at constant composition 
with temperature as parameter. The first three 
cases will deal with gases only, and we will start 
with the relatively simple case represented in 


Fig. 1. This situation is purely hypothetical, and 
probably does not actually exist. According to 
eq. (17), 


— dp 


RT - in( f,/2,f,) 


Therefore will equal z, f, at all pressures up to 
the point where V’, is no longer equal to V, along 
an isotherm. Furthermore, H, will equal H, up 
to nearly the same pressure. This case represents 
a situation where a solution could be ideal for all 
practical purposes over a limited range of con- 
ditions. 


(V,-V,) A 


Schematic diagram. Constant total composition. 


Fie. 1. 
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(Vi-V;) 
p 


Fie. 2. 


Schematic diagram. Constant total composition 

In the situation represented by Fig. 2, no ideal 
solution is possible. Again, this is a hypothetical 
situation which may or may not actually exist. 
It is true that at the temperature T,, (V; — V,) 

0 for a range of pressures down to zero, and 
that f; = z,f, over this pressure range, but H, 
is never equal to H, over this pressure range. 
This may be seen by differentiating eq. (17) with 
respect to temperature after having first divided 
through by T. The differentiation is carried out 
at constant pressure and composition, and gives : 


p 


(V,—V,) dp 


V,) dp- (T) 


p 
-H,= 


oT 


While the first term on the right-hand side of 
this equation is zero over a wide range of pressures 
at T,, the second term will always have a value. 
Thus such a solution may act as an ideal solution 


at this temperature in certain respect’s, but not in 
all respects, and therefore cannot be called an 
ideal solution. 


p 


Fic. 3. Constant total composition. 


Schematic diagram. 


The situation represented in Fig. 3 is commonly 
encountered with gaseous solutions at tempera- 
tures above the critical. Experimental evidence 
shows that some if not all gas mixtures very nearly 
obey Amacat’s law at high pressures (usually 
above several hundred atmospheres), even though 
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at lower pressures, indeed at zero pressure, this 
In such a case the Lewis and 
valid for any 


lay may fail. 
Randall rule is never exactly 
temperature shown, for the integral term of 
eq. (17) always has a value even in the region 
where VF, V,. Moreover, H, does not become 
equal to H, in this region because the value of 
the integral is different for different temperatures. 
However, the ratio, f,/z,f,, becomes constant at 
high pressures for a given temperature, and the 
difference between H, and H, becomes constant 
for the same conditions. This is apparent from 
eqs. (17) and (27). 

Of the three general cases considered above for 
gases, only the first represents a behaiour suitable 
for the exact application of the equations for an 
ideal solution over any finite range. Unfortu- 
nately, this is strictly a hypothetical behaviour, 
not known to be exhibited by any actual gas 
mixture. This analysis lends support to the view 
expressed by Repiicn and Kwone [3] that no 
gaseous mixture can behave as an ideal solution 
unless it is also an ideal gas. This does not mean 
that the use of the equations for an idea! solution 
cannot represent any improvement over the use 
of the equations for an ideal gas, for while 
neither may be correct as applied to actual gases, 
the former certainly may be much more nearly 
so than the latter. 


Constant total composition 


Fic. 4. Schematic diagram. 


The general case represented by Fig. 4 is 
particularly important, for it includes the situation 
which must exist for any ideal liquid solution. If 
one starts with a vapour mixture at the very low 
pressure, p*, constant 
temperature below the critical temperature, the 
dew point will eventually be reached. Further 
compression results in continued liquefaction until 
the bubble point is reached. While it is perhaps 
possible for V, to equal V, in the single-phase 


and compresses at a 


regions, this cannot be the case in the two-phase 
region. 

Consider now a plot of V, the total volume, 
vs. z, for a two-component mixture consisting of 
a and b, with a being the more volatile, at a 
temperature and pressure such that pure 6 is a 
liquid and pure a is a vapour. Such a plot is 
shown in Fig. 5. It should be emphasized again 


ABOVE THE DEW 
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TRE BUBBLE POINT 
POINT 
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AND BUBBLE POINTS 
| 
“a 


mole fraction. Constant temperature 


and pressure. 


Fic. 5. Volume vs. 


that V and z refer to the entire mixture even in 
the two-phase region. From this diagram, using 
the method of tangent intercepts [2], we can 
the 


have. 


relative values that and FV, 


Using this information, we can 


determine 
must 
prepare a graph of total volume, pure-component 
volumes, and partial molal volumes as a function 
binary mixture at 


total 


of pressure for this same 


constant temperature and composition. 


This is shown in Fig. 6. 
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BP 
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Fic. 6. Mixture volume, pure component volumes, and 
partial molal volumes vs. pressure. Constant temperature 
and total composition. 
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Fic. 7. Schematic diagram. Temperature and total com- 
position constant. 


From the values given by Fig. 6, graphs of 
(V, — V,) and (V, — V,) vs. pressure can be 
drawn. These are shown in Fig. 7. Temperature 
and total composition are constant. The curves 
shown are of the shape obtained when Raou.t’s 
law is assumed for the vapour-liquid equilibrium 
relation. In addition additivity of volumes in the 
single-phase regions has been assumed. Thus, we 
have represented here the simplest possible case. 
Clearly it is not possible for (Vv; — V,) to be zero 
at all pressures from zero to the pressure in 
question for a liquid mixture. However, it is 
possible for the plus areas to be equal to the 
minus areas in Figs. 4 or 7, and therefore for the 
integral term of eq. (17) to be zero in the liquid 
region. If this is the case, the liquid may well 
be an ideal solution. The only additional require- 
ment is that this be true not only at the tem- 
perature in question but over a range of tem- 
peratures including the temperature in question ; 
otherwise the heat of mixing would not be zero 
as can be seen from eq. (27). 

It is also possible, of course, for the plus area 
to be greater or less than the minus area, and 
still have a liquid solution for which V; = V,. 
Such a solution could not be ideal, since the 
integral term in eq. (17) would not be zero, and 
therefore the Lewis and Randall rule would 
not be true. Furthermore, in such a case one 
could hardly expect the value of this integral 
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to be the same over a range of temperatures 
even though the liquid volumes remained 
additive, and therefore the heats of mixing could 
not be zero. Such a situation is given a generalized 
representation in Fig. 8. 


Fic. 8. Schematic diagram. Total composition constant. 


There are undoubtedly examples of actual 
liquid solutions for which the volume change of 
mixing is very small over a wide temperature 
range, but which are by no means ideal from the 
standpoint of the Lewis and Randall rule or 
of zero or small beats of mixing. The existence 
of such solutions is easily explained on the basis 
that their vapour-liquid equilibrium relation is 
such as to give a finite value to the integral term 
of eq. (17) through the two-phase region. As a 
matter of fact, the volume changes which occur 
from the mixing of liquids are almost always small 
relative to the volume changes which result when 
mixtures in the two-phase region are formed. Thus, 
it is the nature of the vapour-liquid equilibrium 
relationship which is of major importance in 
determining the magnitude of the heat of mixing 
and the deviation of the activity coefficients 
(f,/2,f,) from unity for liquid solutions. 

This raises the question of what the vapour- 
liquid equilibrium relationship must be in order 
that a liquid be an ideal solution. No completely 
general answer to this question seems to be 
possible; however, if the vapour phase is 
assumed to be an ideal gas, the required relation 
is very nearly Raovu t's law, and if one assumes 
pressure to be without effect on the vapour 
pressures of the pure components, it is exactly 
Raou.t’s law. Clearly, there are many vapour- 
liquid equilibrium relationships which could lead 
to the situation shown in Fig. 7, where the plus 
and minus areas are equal. The conclusion which 
must be drawn from this observation is that 
Raovutt’s law cannot be derived from the simple 
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assumptions that the liquid volumes are additive 
and that the vapour is an ideal gas. The necessary 
assumptions are, that the integral term of eq. (17) 
be zero, and that the vapour be an ideal gas. 
We have been considering the application of 
eq. (17) to the entire system as we have defined 
it. This equation reduces to the Lewis and 
Randall rule for a given component eq. (5) only 
when the integral term is zero. This cannot be 
the case in the two-phase region as can be seen 
by examining Fig. 7. Thus, one cannot apply the 
Lewis and Randall rule to the system as a 
whole if the system consists of two phases. One 
can, however, apply eq. (17), and thus the Lewis 
and Randall rule where it is warranted to the 
individual phases. This often requires the use of 
fictitious states for the pure components, for as 
can be seen from Fig. 7, the integral term in 
eq. (17) cannot be zero in a limited region outside 
the two-phase region when only actual states are 
employed. To avoid this difficulty outside the 


two-phase region, it is customary to use values 
of f, in the Lewis and Randall rule for the pure 


components in the state (liquid or vapour) in 
which the solution exists, even though the pure 
component does not actually exist in this state. 
In such cases the Lewis and Randall rule is 
written f,;=2,f, or f; = depending on 
whether the phase being considered is liquid or 
vapour. In the equations derived in this paper. 
the use of fictitious states has been avoided, 
However, this is not to say that there is anything 
objectionable about their usual and proper use. 

The application of eq. (21) to liquid solutions is 
also of interest. The latent-heat terms are, of 
course, large. For an ideal solution AH is zero. 
and even for non-ideal solutions AH is often 
relatively small. Thus the integral term must 
have a large value. This comes largely from the 
two-phase region, which must be passed in 
reaching the liquid state from the zero-pressure 
state. Thus, AH depends on the difference 
between two relatively large numbers. Represent- 
ative values of AH and AV for liquids at about 
20°C as given by Scumipr [4] are listed in Table 1. 
In each case values are for an equimolal mixture. 
The units are cal/gm of solution for AH and 
ee/gm of solution for AV. 


Van NEss 


Table 1. Representative values of AH and AV. 


System AH AV 
Acetone-Carbon disulphide + 5-00 + 0-0140 
Acetone-Chloroform ~— 5-56 — 0-0020 
Benzene-Methyl alcohol + 1-98 — 0-0014 
Benzene-Carbon disulphide + 1-90 + 0-0005 
Benzene-Ethy] ether 0-00 — 0-0060 
Benzene-Carbon tetrachloride + O25 + 0-0051 
Benzene-Propyl] alcohol + 2-93 0-0052 
Carbon disulphide- 

Carbon tetrachloride + 1-59 0-0053 


These data show clearly that there is no direct 
relationship between AH and AV. However, 
eq. (21) is in accord with such results as is eq. 
(27). Whether or not eq. (21) could be used to 
advantage for the calculation of heats of mixing 
of liquids is questionable. Certainly, very 
accurate volume measurements would be required 
not only in the vapour region but in the two-phase 
region as well. Of course, such measurements are 
being made in the vapour region, but one has the 
added difficulty of assuring phase equilibrium if 
they are extended into the two-phase region. On 
the other hand, the direct determination of 
accurate heats of mixing is also quite difficult. 

We may conclude from this discussion that one 
must be careful not to assume a liquid solution 
to be ideal or nearly so in all respects just because 
it shows negligible or small volume changes of 
mixing in all concentrations, even over a wide 
range of temperatures, for this can be true even 
though the Lewis and Randall rule is not valid. 
and the heat of mixing not zero. To test a liquid 
solution for approach to ideality, one should 
measure not only volume changes of mixing over 
a range of concentration and temperature, but 
also heats of mixing (or temperature changes of 
mixing). 

This paper has dealt only with purely thermo- 
dynamic relationships, and has not been con- 
cerned with the question of ideality from the 
standpoint of molecular structure and inter- 
molecular forces. While it should be possible to 
determine the approach to ideality of a solution 
on its own merits in this way, present knowledge 
does not in general permit it. Thus one must 
have recourse to thermodynamic relationships 
such as are presented here. 
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(1) 


NOTATION 


molal free energy of pure component i at the 
temperature and pressure of the solution 

partial molal free energy of component i in solu- 
tion ; same as the chemical potential, ,;. 

free energy change of mixing at constant tem- 
perature and pressure to form 1 mole of solution 
free energy change of mixing over and above that 
for the mixing of ideal gases 

fugacity of pure component i at the temperature 
and pressure of the solution 


= fugacity of component i in solution 


molal enthalpy of the solution 


= molal enthalpy of pure component i at the tem- 


p= 


perature and pressure of the solution 

partial molal enthalpy of component i in solution 
enthalpy change of mixing at constant tempera- 
ture and pressure to form 1 mole of solution 
(integral heat of mixing) 

molal latent heats of pure a and pure b at the 
temperature of the solution 

pressure 


REFERENCES 


= universal gas constant 
= molal entropy of pure component i at the tem- 


perature and pressure of the solution 


= partial molal entropy of component i in solution 


entropy change of mixing at constant tempera- 
ture and pressure to form 1 mole of solution 
entropy change of mixing over and above that for 
the mixing of ideal gases 

absolute temperature 

molal volume of the solution 

molal volume of pure component i at the tem- 
perature and pressure of the solution 

partial molal volume of component i in solution 
volume change of mixing at constant temperature 
and pressure to form 1 mole of solution 

mole fraction of component i in a liquid phase 
mole fraction of component i in a vapour phase 
mole fraction of component i in the total system, 
whether vapour, liquid, or a combination of both 


= a superscript, designating the zero-pressure state 
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The mechanism of mass transfer of solutes 
across liquid—tiquid interfaces 


Tue two recent papers by J. B. Lewis [2, 3] report data 
on mass transfer across the flat interface between two 
immiscible liquids in a stirred vessel. There appear to 
be several reasons to question the conclusion that mass 
transfer in such a system “ takes place by eddy and not 
molecular diffusion,’’ a conclusion which is at variance 
with current theory. 

The common picture of the mechanism of mass transfer 
to a phase boundary is one of transport by eddy motion 
in the flowing or agitated bulk of the fluid, with eddy 
motion in a direction normal to the surface damped out 
near the surface. Lack of perpendicular velocity com- 
ponents, resulting from this damping, suggest that in the 
limit, as the surface is approached, molecular motion 
is the only possible mechanism of mass transfer. It is 
difficult to see how lateral motion of the surface, as in the 
ease of two stirred immiscible liquids, would change this 
picture. 

Molecular diffusion is notoriously slow as compared 
with eddy diffusion, and the effect of the molecular 
diffusion coefficient D is not significant in the stirred bulk 
of the fluid. As the surface is approached the rate of 
transfer should, in the limit, be proportional to D. It may 
be expected, then, that the over-all process will be pro- 
portional to D raised to some power between zero and 
unity. Empirical correlations of mass transfer data 
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frequently involve D”", where n is in the vicinity of 0-5. 
Variations in turbulence affect both eddy transport and 
the effective thickness of the region near the surface where 
D is of importance, so GrLLILanp [1], for example, found 
no evident variation of n as the stream velocity was 
varied several fold. The Danckwerts surface-renewal 
model, likewise, calls for a constant value of n of 0-5. 

Because of the large number of variables involved, 
Lewis’ task of correlating his extensive data was admittedly 
difficult, but it was not made easier by the fact that his 
results were not reproducible. 

In a number of cases the values of the mass transfer 
coefficient & differed by more than two-fold in presumably 
identical tests. In spite of the fact that smoothed data 
were used, the final correlation (Fig. 4 of the first paper) 
shows a wide spatter of points, and the furfural data lie 
more than 100 per cent above the best line. The resulting 
correlation is curious in that the viscosity of the “ other 
phase” cancels. Over a wide range the line of Fig. 4 
has a slope of nearly unity so that the viscosity of the 
phase under consideration also cancels. The result, 
roughly, is a correlation of k with stirrer speeds. 

Under the circumstances, it is not unreasonable to 
believe that other methods of correlation might serve 
equally well, and that a correlation involving D might 
have been obtained. Too few reliable values of D are 
available to attempt this properly, but the possibility 
can de demonstrated. 
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Replotting the original data in the manner of Lewis, 
Fig. 3, values of k were interpolated for Reynolds numbers 
in both phases of 2,000 and 4,000. This smoothing proce- 
dure is quite uncertain, but a few values good to perhaps 
80 per cent were obtained for those systems for which 
Lewis reports [4] values of D. These are shown on the 
accompanying logarithmic graphs (Figs. 1 and 2) ask vs. D, 
with lines of slope 0-5 sketched in. This fragmentary 
correlation does not, at least, do violence to the idea that 
k depends in some way on D, and it would not be unreason- 
able to conclude that & is proportional to the square root 
of D, as found by other workers. The coefficient k appears 
to increase regularly with the viscosity of the organic phase, 
though the values of k shown are for diffusion in both 
organic and aqueous phase. Perhaps the more viscous 
organic liquids caused increased shear at the interface, 
with a resulting increase in k in each phase, quite apart 


from the effect of Reynolds number. It cannot be claimed 
that this approach to a correlation is comprehensive or 
even valid-only that one involving D cannot be ruled 
out, and that Lewis’ conclusion that molecular diffusion 
plays no part must be questioned. 

If the correlation presented in the first paper is suspect, 
then the analysis of the data of the second paper dealing 
with mass transfer of a third component, is even more 
Quite apart from the possibility of other 
correlation that the 
reproducibility of the data and the deviations of experi 
mental k’s from Lewis’ general correlation would throw 


questionable. 


procedures, it would seem poor 


serious doubt on the reliability of the calculated values 


of k,,, and perhaps even of the existence of a “ slow 
interfacial reaction.” 
Tuomas K. Suerwoop 


Massachusetts Institute of Technology, Cambridge, Mass. 
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In reply to Professor Sherwood 


ALTnouGn the presence of a laminar boundary layer at 
solid-liquid interfaces has been conclusively demonstrated 
experimental difficulties hitherto prevented the 
direct investigation of flow patterns at gas-liquid and 
liquid-liquid interfaces. Nevertheless it 
postulated, from analogy with the solid-liquid surface, 


have 
is commonly 


that similar boundary layers are always obtained, under 
all stirring conditions, at inte -faces This 
postulate has been accepted as correct because it is possible 
to correlate experimental data on the basis of a molecular 
It is 
not at all clear why this postulate, which assumes zero 


these also. 


diffusional process being the rate controlliny step. 


perpendicular velocity at gas-liquid and_ liquid-liquid 


interfaces, arose since in these cases there is no rigid 
immobile surface as at the solid-liquid bondary. Again 
it has recently been shown by Danckwerrts [1], in 
connection with gas absorption, that different transport 
mechanisms will nevertheless all lead to transfer equations 
of essentially the same form. Although Kisuinevski [2)} 
has criticised DaNcKwerts’ own postulated mechanism 
and has shown that absorption could well be entirely 
turbulence controlled he again derives a transfer equation 
of the customary type. Thus it is clear that ability to 
correlate experimental data according to an equation 
derived from a diffusional layer theory is in itself no proof 
of the validity of the theory. 

In my work the most direct evidence for the absence 
of laminar boundary layers, under the particular stirring 
conditions used, was that the transfer coefficients were 
independent of the relative directions of rotation of the 
two stirrers. Most experiments were carried out with 


contra-rotating stirrers, to produce maximum shear near 
the interface, which was then always disturbed by small 
ripples. In 10 experiments both stirrers were rotated in 
the same direction giving a very smooth interface with 
minimum shear but there was no significant difference in 
the transfer coefficients. It was possible to start an experi- 
ment with contra-rotating stirrers and to reverse one 
stirrer midway through the experiment without a break 
appearing in the (linear) plot of log (c,, — c) versus time 
Clearly if laminar layers had been present they could not 
have been of the same thickness in the two casei. On the 
other hand if transfer was by eddy diffusion the important 
factor would be the component of velocity perpendicular 
to the interface, which would be independent of the 
direction of stirring. 

A further check on the validity of a diffusion il process 
was made by measuring transfer rates at different tempera- 
tures. At 50° it was clear, in spite of the experimental 
scatter, that there was not the large increase in transfer 
rate to be expected if the rates were diffusion controlled 
Again the much greater retardation produced by rigid 
interfacial films on the transfer rate, as compared with 
mobile films, can be explained by the assumption that in the 
former case the rigid interface allows laminar layers to be 
built up whereas the mobile film has no effect. 

The scatter of results in Fig. 4 arises mainly from the 
somewhat different coefficients observed at similar values 
of the group (Re, + Rez n2/m,) which were however 
obtained from different combinations of the two Reynolds 
numbers. This scatter is much les; when the two Reynolds 
numbers are equal. An idea of the average experimental 
scatter can be obtained from the data plotted in Fig. 3. 
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It is clear that there is at least one further term missing The criticism that over a wide range of Fig 4, the two 
from the correlation but provided the two Reynolds viscosities both cancel is hardly valid since this merely 
numbers are approximately equal then the data can be demonstrates once again that the process is turbulence 
predicted to sufficient accuracy by the existing graph. controlled. 

The experiments reported in Part II were carried out as Chemical Engineering Division A.E.R.E., 
far as possible with the two Reynolds numbers equal. Harwell, Didcot, Berks. 


J. B. Lewis 
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ANNOUNCEMENT 


International Congress on Catalysis 


HONORARY CHAIRMEN 
Mea. Evcene J. Houpry, Sim Eric Ripear, Pror. H. Tayior 


Philadelphia, Pa., September 10-14, 1956 


Four major symposia are envisaged, each of which will deal with a variety of subjects. 


1. Chemistry and Physics of Solid Catalysts (Catalyst Structure, Electron Transfer Processes in Solids, 
Poisoning, Reaction Kinetics). 


Homogeneous Catalysis and Related effects (Catalysis by Radicals, Catalysis by Ions, Biocatalysis, 
Catalysis and Radiation, Photosynthesis). 


Surface Chemistry and its Relation to Catalysis (Corrosion, Flotation, Adsorption, Electrode Processes, 
Colloids). 


Techniques and Technology of Catalysis (Laboratory Techniques for Studying Reactions, Catalysis 
Preparation, Catalytic Reactions of Technological Importance). 


All papers on the programme will be circulated in the form of preprints, and brief summaries only will be presented at 
the meetings, to assure maximum opportunity and time for discussion. A few major topics will be presented in more 
detailed lectures by invited speakers. English, French and German are the official languages of the Congress. 


Prospective contributors are requested to submit the titles and a two-hundred-word abstract for their papers not later 
than January 1, 1956. The full paper must be ready for preprinting by March 1, 1956. All papers submitted for 
presentation at the Congress will be reviewed by a committee before acceptance. 


Correspondence concerning the programme should be sent to Dr. A. Farkas, Chairman, Programme Committee, 
International Congress on Catalysis, Barrett Division, Allied Chemical & Dye Corporation, S. Ridgeway Avenue, 
Glenolden, Pennsylvania. All other communications should be addressed to Dr. H. Heinemann, Executive Secretary, 
International Congress on Catalysis, ¢/o Houdry Process Corporation, P.O. Box 427, Marcus Hook, Pennsylvania. 
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